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ABSTRACT 

The  omni-present  low  frequency  wave  motion  (30-300  sec)  contains  a  substantial 
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the  shoreline).  The  modulation  is  identified  in  the  analysis  of  field  data  as  a  positive 
correlation  between  the  long  waves  and  the  wind  wave  envelope  near  the  shoreline. 
Considering  oblique  incident  waves,  a  steady  longshore  current  showing  a 
non-vanishing  current  at  the  shoreline  is  found  as  a  result  of  the  non-linear  interaction 
between  monochromatic  incident  and  infragravity  waves.  An  analytical  solution 
describing  the  unsteadiness  of  the  longshore  current  is  developed.  Applying  the  derived 
longshore  current  solution,  longshore  sediment  transport  is  reformulated  to  include  the 
infragravity  waves,  giving  improved  comparisons  with  field  measurements. 
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I.  INTRODUCTION 

Sea  waves  are  generated  by  wind,  and  as  they  propagate  towards  the  shore  they 
deform  due  to  shoaling  and  may  refract  if  they  are  oblique  to  the  bottom  contours  or  if 
they  meet  a  current.  The  waves  grow  in  height  until  they  become  unstable  and  break, 
dissipating  their  energy.  The  processes  of  breaking  waves  on  a  sloping  bottom  make 
the  surf  zone  an  extremely  dynamic  area  where  different  observable  phenomenac  take 
place,  such  as  longshore  current,  sediment  transport  and  wave  runup  (setup  and 
swash).  Due  to  this  complexity,  the  proper  modeling  of  wave  motion  in  the  ncarshore 
zone  has  been  the  goal  of  many  investigators,  and  yet  much  remains  to  be  done  before 
a  satisfactory  understanding  is  achieved. 

The  surf  zone  can  be  characterized  dynamically  by  three  regions:  the  outer.thc 
inner,  and  the  runup  region  (Svendson  et  al.,  1978),  (see  Fig  1.1).  The  outer  region 
follows  the  breaking  point  and  is  dominated  by  incident  breaking  waves.  The  breaking 
waves  are  of  either  the  plunging  or  spilling  type,  and  are  characterized  by  a  rapid 
transition  of  wave  shape  with  a  horizontal  surface  roller.  In  the  inner  region,  the  wave 
forms  break  down  into  smaller  scales  of  random  and  turbulent  nature  and  eventually 
become  similar  to  moving  bores,  or  hydraulic  jumps.  The  run-up  region  is 
characterized  by  wave  setup  and  swash  oscillation  running  up  and  down  the  beach 
face.  In  all  regions,  the  wave  height  and  the  water  motion  are  strongly  locally 
controlled  by  the  depth.   The  depth  controlled  waves  arc  referred  to  as  "saturated." 

In  general,  surf  zone  dynamics  are  found  to  differ  dramatically  between  rellectivc 
(steep)  beaches  and  dissipative  (gentle  slope)  beaches  (Bradshaw,  1980).  Typical 
features  of  the  reflective  beaches  are  high  reflection  of  incident  waves,  collapsing  or 
surging  breakers  on  the  beach  face,  and  subharmonic  resonances  (e.g.  Guza  and  Davis, 
1974;  Guza  and  Inman,  1975;  Huntley  and  Bowen,  1975b,  1978;  Wright  ct  al.,  1979). 
Dissipative  beaches  arc  dominated  by  spilling  breakers,  pronounced  low  frequency 
oscillations  and  the  presence  of  multiple  parallel  longshore  bars  (e.g.  Huntley  and 
Bowen,  1975a;  Short,  1975;  Sasaki  and  llorikawa,  1975;  Ilolman  ct  al.,  1978;  Symonds 
ct  al.,  1982).  This  study  will  focus  on  the  dynamics  of  dissipative  beaches  in  which  the 
inner  region  is  assumed  to  be  a  transition  region  between  regions  dominated  by 
incident  breaking  waves  and  low  frequency  swash. 
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Figure  1.1     Surf  zone  regions  and  coordinate  system. 
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In  the  past  decade,  considerable  attention  in  the  literature  of  bcacli  and 
nearshore  dynamics  has  been  devoted  to  the  kinematics  of  low  frequency  waves.  These 
waves,  which  have  periods  ranging  from  30-300  sec,  are  referred  to  as  "infragravity 
waves."  The  observations  show  that  the  dominant  wave  motion  closer  to  the  shore  is 
not  normally  at  the  incident  wave  frequency  band.  It  is  also  found  that  the  low 
frequency  wave  energy  in  the  very  nearshore  region  contains  a  substantial  fraction  of 
the  total  wave  energy  (Inman,  1968  a,b;  Suhayda,  1974;  Goda,  1975;  Huntley  and 
Bowen,  1975b;  Sasaki  and  Ilorikawa,  1975,1978;  Huntley,  1976;  Wright  et  al., 
1978,1979,1982;  Bradshaw  19S0;  Holman,  1981;  Huntley  ct  al.,  1981;  Guza  and 
Thornton,  1982, 1985a,  1985b;  and  others),  (see  Fig  1.2  for  example).  The  valley  in  the 
energy  spectra,  where  the  energy  level  reaches  minimum,  suggests  an  empirical  cutoff 
frequency  around  0.03  Hz  that  separates  the  two  wave  bands.  Therefore,  a  complete 
description  of  the  wave  field  in  the  surf  zone  should  include  two  relatively  separated 
bands  of  wave  frequencies  :  (1)  incident  sea-swell  waves  (short  period  waves)  with 
periods  of  1-30  sec,  and  (2)  infragravity  waves  (long  period  waves)  with  periods  of 
30-300  sec,  which  no  longer  can  be  neglected  inside  the  surf  zone. 

Prior  to  the  1950's,  the  incident  wave  characteristics  were  the  main  concern  when 
measuring  and  analyzing  waves  in  the  nearshore  zone.  Munk  (1949)  and  Tucker  (1950) 
were  the  first  to  point  out  the  existence  and  significance  of  low  frequency  motion  and 
called  the  phenomenon  "surf  beat."  More  recently,  with  renewed  interest  in  nearshore 
processes  and  with  major  improvements  in  the  instrumentation,  it  has  become  apparent 
that  the  low  frequency  wave  energy  can  be  as  important  as  the  incident  wave  energy 
when  approaching  the  shoreline.  Another  motivation  for  the  study  of  the  infragravity 
waves  is  the  coupling  between  the  infragravity  length  scales  (order  of  hundreds  of 
meters)  and  the  length  scales  associated  with  most  of  the  rhythmic  beach  features  such 
as  sand  bars,  beach  cusps  and  crcscentic  bars.  Many  investigators  have  related  the 
beach  morphology  changes  to  the  existence  of  low  frequency  motions  by  which  a 
mechanism  of  forming  and  maintaining  such  topographic  features  is  provided  (Bowen 
and  Inman,  1971;  Guza  and  Inman,  1975;  Bowen,  1980;  Holman  and  Bowen,  1982; 
Bowen  and  Huntley,  1984;  and  others). 

In  analyzing  surf  zone  data,  there  are  often  anomalous  patterns  of  water  motion 
closer  to  the  shoreline  when  compared  with  present  short  wave  models,  for  example, 
the  spatial  distribution  of  wave  height  measured  in  both  the  laboratory  and  field  often 
shows  unexpected  wave  height  at  the  shoreline  (e.g.    Suhayda  and  Pettigrcw,   1977: 
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Figure  1.2     Cross-shore  velocity  spectra  on  20  Nov.    1978,  Torrcy  Pines 
(alter  Gu/.a  and  Thornton,  1985b). 
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Battjes  and  Stive,  1985),  (sec  for  example  Fig  1.3).  Bowen,  ct  al.  (1968),  in  measuring 
the  setup,  found  a  residual  wave  height  always  present  at  the  shoreline  and  modified 
their  setup  formula  at  the  beach  face  to  lit  the  data.  Suhayda  and  Pcttigrcw  (1977) 
measured  the  average  wave  velocity  across  the  nearshore  zone  and  found  it  to  decrease 
shoreward  until  the  sea-swash  limit  of  the  run-up  zone  was  reached,  and  then  to 
increase  onshore,  supporting  Waddell's  measurement  (1973).  Thus,  the  wave  height 
and  celerity  do  not  vanish  at  the  vicinity  of  the  shoreline,  invalidating  the  widely 
applied  short  wave  saturation  assumption. 

The  assumptions  of  wave  saturation  and  steady  state  conditions  for  both 
monochromatic  and  random  waves,  results  in  a  zero  current  at  the  shoreline.  In 
addition,  the  steady  longshore  current  assumption  is  often  unrealistic,  especially  with 
the  presence  of  irregular  oscillating  waves  (Wood  and  Meadows,  1975).  This  is  often 
seen  in  field  data  that  show  longshore  velocity  oscillations  at  longer  wave  periods  than 
for  incident  waves  (e.g.  Holman  and  Bowen,  1984;  Guza  and  Thornton.  1985a; 
Oltman-Shay  and  Guza,  1986).  Using  an  empirical  orthogonal  function  (EOF)  analysis 
on  64  data  sets  covering  a  wide  variety  of  wave  conditions,  Guza  et  al.  (1986)  obtained 
objectively  best  fit  to  the  longshore  current  distribution.  The  analysis  showed  the 
"classical  parabolic"  shape  for  the  longshore  current  distribution  (Bowen,  1969)  over 
most  of  the  surf  zone,  but  with  significant  residual  velocity  at  the  shoreline  (Fig.  1.4). 
In  the  analysis  of  the  same  data,  Wu  et  al.,  (1985)  found  a  poor  agreement  between 
longshore  current  theory  with  velocity  going  to  zero  at  the  shoreline  and  data  closer 
to  the  shoreline.  All  present  models  predict  zero  current  at  the  shoreline,  and  no  model 
is  available  to  explain  the  unsteadiness  of  longshore  current. 

Improving  the  longshore  current  description  is  important  because  it  provides  the 
driving  means  of  transporting  sediment  alongshore  (Longuct-IIiggins,  1972).  The 
available  longshore  sediment  transport  models  (e.g.  Bijker,  1971;  Thornton,  1973; 
Komar,  1977;  Sawaragi  and  Dcguchi,  1979;  Bailard,  1981),  all  show  a  longshore 
transport  maximum  closer  to  the  brcakerlinc,  which  agrees  with  the  longshore  current 
distribution.  The  shortcoming  of  these  models  is  their  failure  to  predict  any  moving 
sand  at  the  shoreline  regardless  of  the  wave  conditions.  Zenkovich  (1960)  measured 
the  suspended  sediment  transport  distribution  across  the  surf  zone  using  lluorcsccnt 
sand  tracers.  Sand  transport  was  found  greater  over  bars  where  energy  dissipation  is  a 
maximum  due  to  wave  breaking  and  also  at  the  shoreline,  (see  Fig  1.5).  Measuring 
longshore  sand  transport  distribution  using  sand  tracers,  Kraus  et  al.  (1981)  found  a 
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Figure  1.3     Root  mean  square  wave  height  HTrm  vs  o  fish  ore  distance  x. 

Solid  curve  represents  a  short  wave  modcrand  ffiU  stars  arc  laboratory  data 

showing  non-zero  wave  height  at  the  shoreline  (alter  liattjes  and  Stive,  1985) 
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distinct  bimodal  distribution  having  maxima  at  the  brcakerline  and  near  the  shoreline. 
The  same  results  were  obtained  by  Downing  (1984)  when  measuring  the  suspended 
longshore  sediment  distribution.  Abdclrahman  (1983)  compared  longshore  sediment 
transport  at  Santa  Barbara  beach  with  his  predicted  longshore  sediment  transport 
model  based  on  the  widely  accepted  energetic  concept  by  Bagnold  (1963).  The  model 
did  not  predict  the  large  transport  in  the  vicinity  of  the  shoreline  since  theoretical  short 
wave  description  gives  zero  amplitude  at  the  shoreline.  It  was  suggested  a  primary 
cause  of  these  difTerences  was  the  presence  of  infragravity  waves,  which  have  not  been 
taken  into  consideration  before. 

In  recognition  of  the  importance  of  the  infragravity  waves  in  the  surf  zone,  an 
obvious  next  improvement  in  the  description  of  waves  kinematics  is  to  combine  the 
infragravity  waves  with  the  incident  sea-swell.  As  a  result,  the  incident  waves  would  be 
modulated  in  amplitude,  wavenumber  and  direction.  The  purpose  of  this  dissertation  is 
to  improve  the  understanding  of  nearshore  processes  by  generating  and  applying  such  a 
composite  wave  description. 

As  will  be  seen,  the  interaction  of  the  short  waves  with  the  long  waves  results  in 
energy  not  only  at  their  respective  frequencies,  but  energy  transferred  to  side  bands  of 
the  incident  waves  at  the  sum  and  difference  frequencies.  It  is  the  energy  and 
momentum  at  the  side  bands  that  is  primarily  responsible  for  the  inner  surf  zone 
dynamics.  The  dynamical  changes  are  demonstrated  by  the  successful  treatment  of 
using  the  first  order  wave  description  in  the  momentum  flux  equation  to  obtain  a 
second  order  dynamical  effect.  However,  inclusion  of  infragravity  waves  does  not 
mean  introducing  more  independent  parameters  since  the  infragravity  wave  energy  is 
assumed  related  to  the  incident  wave  field. 

In  the  following  chapter,  a  description  for  the  short  and  long  period  waves  is 
presented.  The  short  period  waves  arc  described  by  linear  (Airy)  theory,  while  the  long 
period  waves  on  a  sloping  bottom  are  described  by  either  cross-shore  standing  waves 
(2-D)  or  edge  waves  (3-D).  In  chapter  3,  the  changes  in  the  amplitude  and 
wavenumber  of  the  short  waves  superposed  on  long  waves  are  discussed.  An  analytical 
expression  describing  the  modulated  surface  elevation  of  the  short  waves  is  developed 
alon°  with  their  encrsy  and  momentum.  Changes  in  the  wave  kinematics  are 
considered  with  the  long  period  waves  treated  as  a  slowly  varying  current.  In  chapter 
4,  narrow  banded  wave  data  are  analyzed  to  examine  the  hypothesis  of  short  wave 
modulation  by  the  long  waves.    The  dynamical  effects  of  considering  the  infragravity 
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Figure  1.4     Longshore  current  distribution  using  empirical  orthogonal  function 

analysis.   The  first  cienfunction  contains  most  of  the  variance 

(after  Guza  et  al.,  llAS6). 
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waves  inside  the  surf  zone  are  investigated.  Analytical  models  for  calculating  wave 
setup,  longshore  current  and  sediment  transport  are  developed  in  chapter  5.  An 
analytical  solution  to  the  unsteady  longshore  current  is  derived.  Model  validation  and 
comparison  with  other  models  are  included  in  chapter  5.  The  conclusions  of  this  study 
are  presented  in  chapter  6. 
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Figure  1.5     Distribution  of  longshore  velocity  and  sediment  transport 
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II.  WAVES  IN  THE  NEARSHORE  ZONE 

A  wave  record  of  incoming  waves  usually  contains  short  period  waves  (1-30  sec) 
and  long  waves  (30-300  sec).  To  study  first  order  modulations  of  short  waves 
superimposed  on  long  waves,  the  principle  of  superposition  is  applied.  Therefore,  the 
non-linear  wave  theories  are  excluded  in  this  study.  A  chosen  Cartesian  coordinate 
system  is  illustrated  in  Fig.  1.1  where  the  positive  X-axis  is  directed  offshore,  the 
positive  Z-axis  is  vertically  upward  and  Y  is  alongshore.  This  coordinate  system  is 
convenient  for  describing  the  offshore  decay  of  the  standing  long  waves.  Linear 
descriptions  of  both  sea-swell  and  infragravity  waves  are  discussed  separately.  Short 
waves  are  derived  from  linear  wave  theory  over  a  flat  bottom,  then  using  shallow  water 
approximation  to  describe  short  waves  in  the  ncarshorc  zone.  Infragravity  waves  are 
derived  from  shallow  water  equations  on  a  sloping  bottom. 

A.       SHORT  PERIOD  WAVES 

In  general,  waves  in  a  viscous  fluid  propagate  over  irregular  topography  of 
varying  permeability.  In  most  cases,  the  main  body  of  the  fluid  motion  is  reasonably 
assumed  irrotational  since  the  viscous  effects  are  usually  limited  to  a  thin  boundary 
layer  near  the  surface  and  the  bottom.  Surface  waves  arc  considered  a  boundary  value 
problem. 

Linear  wave  theory  (small  amplitude  waves)  is  considered  a  first  order 
approximation  to  the  theoretical  description  of  wave  behaviour.  It  is  assumed  that  the 
wave  amplitude  ,  a,  is  very  small  compared  to  both  the  wavelength  L  and  the  local 
water  depth  h,  (a  <  <  L  and  a<<  h).  The  fluid  is  assumed  incompressible, 
homogenous,  inviscid  and  also  irrotational,  by  which  a  velocity  potential  (p  should  exist 
to  represent  the  field  of  the  flow: 


u      =      -    Vo  (2.1) 
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where  u  is  the  velocity  vector.    liquation  2.1  satisfies  the  continuity  equation  V  .u  =  0, 
by  which  the  governing  differential  equation  (Laplace  equation)  is  obtained: 

(   "h   i  z   1  n 
V2$      =      0  (2.2) 

(      _,r>      <       X      <      °° 


With  the  proper  kinematical  and  dynamical  free  surface  and  rigid  flat  bottom  boundary 
conditions,  the  method  of  separation  of  variables  is  used  to  solve  equation  2.2  Then, 
the  velocity  potential  is  given  by 


_  ag  cosh  k(h+z)  cos(kx  +  ut) 


(jj 


sinh  kh 


where  a  is  the  wave  amplitude,  CO  =  (2lt/T)  is  the  radial  frequency,  and  k  =  (2rt/L)  is 
the  wavenumber.    The  surface  elevation  i|  is  periodic  in  both  space  and  time. 


n      =     acos(k.x.    +ojt)     ,        i  -   1,2 


(2.4) 


where  i  =    1,2  refers  to  conditions  in  the  X  and  Y-dircctions.    The  description  of  the 
wave  motion  is  completed  by  the  dispersion  relationship, 


2 
oj  gk   tanhkh  (2.5) 


which  arises  from  the  dynamical  free  surface  boundary  conditions.    The  horizontal 
velocity  under  the  wave  is  given  by 


cosh    k    h+z  ,,        .       .  » 

aj ■■ r — r-r COS       (kx     +     (jjt)  n  sx 

sinh    kh  (-•») 
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and  the  vertical  velocity 


T7             *      sinh    k(h+z)       .     ,.  ^.  (2  7) 

W      =      a"      sinh    kh sm(kx   +   a)t)  ^'> 


Swell  waves  propagating  from  deep  water  are  nearly  sinusoidal  with  long  low 
crests.  As  waves  approach  shallow  water,  they  start  to  feel  the  bottom  and  deform  at 
relative  depth  h/L  <  1/2.  The  wavelength  decreases  while  the  wave  height  generally 
increases,  and  the  wave  period  remains  invariant.  The  hyperbolic  functions  have 
convenient  deep  and  shallow  water  asymptotes  by  which  an  approximate  solution  can 
be  obtained.  In  shallow  water,  the  linear  gravity  waves  are  non-dispersive  since  the 
dispersion  relationship  reduces  to 


J      =      gk2h  (2.8) 


and  the  horizontal  and  vertical  velocities  simplify  to 


u 


=      a    /g/h    cosfkx    +   cot)  (2.9) 


W      =      aco  (1  +£-)    sin  (kx   +   tot)  (2  10) 


On  dissipativc  beaches,  long  crested  swell  waves  break  when  the  wave  steepness 
increases  and  the  velocity  of  the  water  particles  at  the  crest  exceeds  the  phase  speed  of 
the  wave  form.  In  the  inner  surf /.one,  bores  retain  their  relatively  long  crested  form 
and  progressively  decrease  in  height  as  they  advance  onshore.  Due  to  turbulence  and 
non-linearities  associated  with  the  breaking  process,  the  motion  of  the  water  particles 
after  breaking  can  no  longer  be  described  analytically.    As  a  first  approximation,  the 
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wave  or  bore  height,  II,  at  any  location  inside  the  surf  zone  is  limited  by  the  local 
depth  according  to  the  saturation  assumption,  i.e. 


H      =      2yh 


(2.11) 


where  y  is  a  constant  of  order  unity.  This  relation  holds  for  spilling  breakers,  the  most 
common  breaker  type  on  natural  beaches.  The  breaker  index  y  is  important  in 
describing  the  breaking  waves  used  in  longshore  current  and  sediment  transport 
models.  This  linear  relationship  is  expected  on  dimensional  grounds  (Longuet-I  liggins, 
1972)  and  gives  a  reasonable  description  of  laboratory  results  lor  monochromatic 
waves  (e.g.  Galvin  and  l;agleson,1965)  and  field  measurements  during  conditions  of 
both  broad  and  narrow  spectral  distributions  of  wave  energy  (Thornton  and  Guza, 
1982),  (fig.  2.1). 

Very  close  to  the  shoreline  in  the  run-up  region,  the  fluctuating  component  of  the 
wave  runup  about  the  mean  water  level  is  known  as  "swash  oscillation."  Miche  (1951) 
hypothesized  that  monochromatic  waves  within  the  surf  zone  may  be  composed  of 
both  a  progressive  component,  which  is  dissipated  onshore  and  has  zero  shoreline 
amplitude,  and  a  standing  component,  which  has  its  maximum  amplitude  at  the 
shoreline.  Many  laboratory  studies  (iMoraes,  1970;  Dattjes,  1974;  Guza  and  Bowen, 
1976;  Van  Dorn,  1976)  report  that  the  shoreline  swash  amplitude  a  is  limited  by  the 
swash  parameter  c    , 

2 
a    u) 

es      =      -S-y  (2.12) 

gS 

where  tan  [J  is  the  beach  slope  and  cs  varies  from  one  experiment  to  the  other  in  the 
range  from  1  to  3.  Guza  and  Thornton  (1982)  attributed  this  variation  to  the 
ill-defined  nature  of  the  backwash  (run-down).  Hunt  (1959)  related  the  total  runup  Rv 
,  defined  as  the  sum  of  setup  and  half  the  swash  amplitude,  to  the  surf  similarity  £,q 
with  an  empirical  constant  M: 
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c  2.1     Wave  heights  inside  the  surf  zone  defining  breaker  inde? 
and  showing  saturation  (after  Thornton  and  Guza,  1982). 
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Rv  /Hsig  =  M  40  =  M  tan  P  /  V  (HQ  /  LQ)  (2.13) 


where  H„ia  is  the  significant  wave  height  and  the  subscript  0  denotes  the  deep  water 
conditions. 

Huntley  et  al.  (1977)  showed  a  swash  "saturation"  condition  at  incident  wave 
frequency,  i.e.  with  increasing  wave  height,  the  steady  set-up  will  increase  but  the 
swash  amplitude  will  not.  Guza  and  Thornton  (1982)  found  a  co  spectral  decay  at 
incident  wave  frequencies,  indicating  that  energy  levels  are  independent  of  incident 
wave  height,  although  Huntley  ct  al.  (1977)  have  shown  co  dependency.  The  swash 
saturation  condition  is  found  to  be  true  for  both  monochromatic  and  random  swash 
fluctuations  (Guza  et  al.,  1984).  Swash  at  low  frequency  is  found  to  be  unsaturated, 
i.e.  the  swash  amplitude  increases  with  increasing  incident  wave  height.  Most  of  the 
low  frequency  energy  is  found  to  have  infragravity  wave  periods.  Thus  swash  spectra 
show  a  saturated  region  at  incident  wave  frequency  and  an  unsaturated  region  at  low 
frequencies  that  dominates  the  whole  spectrum. 

To  evaluate  the  amount  of  runup,  swash  studies  have  been  mainly  concerned 
with  measuring  the  swash  amplitude  at  the  shoreline  without  putting  much  emphasis 
on  the  generation  and  types  of  the  dominant  low  frequency  motions.  It  is  of  interest 
now  to  investigate  the  forms  of  this  low  frequency  water  motion  that  give  maximum 
amplitude  at  the  shoreline. 


B.       LONG  PERIOD  WAVES 

The  long  period  waves  presented  in  this  work,  are  limited  to  waves  having  a  time 
scale  of  30-300  seconds  where  gravity  is  the  primary  restoring  force. 

1.  Infragravity  Wave  Models 

The  principal  wave  types  that  may  contribute  to  the  low  frequency  motion  on 

beaches  can  be  summarized  by  their  mathematical  formulation  as  follows: 

(i)       2-D  model  :  the  wave  motion  is  mainly  in  the  cross-shore  direction  and  is 
known  as  surf  beat  or  "leaky  modes." 

(i)       3-D  model  :  the  wave  motion  has  longshore  variation  and  is  known  as  edge 
waves  or  "trapped  modes". 
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The  mathematical  forms  of  both  waves  are  derived  from  the  linear  shallow 
water  wave  equation  on  a  sloping  bottom 


2fl  -  g{f_(h  3»)    +|_(h|t))     =     0  (2-14) 

.  1 2         ^    3  x         ox'  dy         oY 


In  the  2-D  model  ,  the  wave  energy  is  reflected  offshore  and  the  solution 
obtained  is  a  standing  wave  expressed  in  terms  of  the  zero  order  Bessel  function  of  the 
first  kind,  Jq  ,  (Lamb,  1932;  Frcdrichs,  1948) 


$(x,t)      =    -^2-  JQ(X)sin   u)t  (2.15) 

where  "/  =  (4«o  x/  g  tan  P)  '  and  a  is  the  amplitude  at  the  shoreline.  It  is  noted  that 
the  Bessel  function  of  the  second  kind,  Yq  ,  is  excluded  since  it  is  unbounded  as  waves 
approach  the  shoreline.  The  surface  elevation  and  the  associated  orbital  velocities  are 
given  by 


n(x,t)      =      a  JJXjcos  ut 


(2.17) 


u(x,t)       =      a    /g/h   J.  (X)    sinut 


W(x,t)       =      -  au  JQ(X)  sin    cot  (2.18) 


In  the  3-D  models,  the  edge  waves  are  free  waves  that  propagate  along  the 
coast  with  their  energy  trapped  onshore  by  refraction,  i.e.  they  do  not  radiate  energy 
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offshore.  Stokes  (1846)  was  the  first  to  provide  an  analytical  solution  to  describe  edge 
waves  using  a  small  amplitude  wave  assumption.  Eckart  (1951)  solved  equation  2.14 
analytically  in  terms  of  Laguerre  polynomials  L  (2k  x).  The  solution  of  progressive 
cdec  waves  is 


-k    x 
4,     =     23-  e     y     Ln(2  kxjcoslky  -   cut)  (2  ,9) 


n 


where  ky  is  the  longshore  wavenumber,  and  n  is  an  integer  number  which  represents 
the  edge  wave  mode  number.  Each  discrete  mode  contains  finite  energy,  and  the  n-th 
mode  edge  wave  will  have  n  zero-crossings  after  which  the  amplitude  decays 
exponentially.  Edge  waves  obtained  by  Eckart  (1951)  must  satisfy  the  dispersion 
relationship, 


g  k     (2n+l)    tan    0  (2.20) 


which  requires  (2n  +  l)tanJ3  <  <   1  to  insure  that  the  solution  achieves  its  limiting  value 
as  k      x  ->-  oo  while  still  in  shallow  water. 

Ursell  (1952)  was  able  to  obtain  a  set  of  "exact"  edge  wave  solutions  on  a 
sloping  bottom  using  the  small  amplitude  wave  theory  without  recourse  to  the  shallow 
water  approximation.    His  solution  satisfies  a  slightly  different  dispersion  relationship 


u*      =      g  k    sin(2n+l)B  (2.21) 


and  the  solution  lies  within  the  range  ((2n+  1)|3  )  <  ju/2  to  insure  a  cutoff  mode  for  a 
given  beach  slope.  Ursell  (1952)  showed  that  Stokes'(1846)  edge  wave  solution  is  only 
the  zero  mode  of  his  solution.  For  gentle  beach  slopes  and  low  mode  numbers,  both 
dispersion  relationships  2.20  and  2.21  arc  in  agreement,  suggesting  that  the  shallow 
water  approximation  is  appropriate  under  these  conditions,  (Guza  and  Davis,  1974).    It 
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is  noted  that  the  approximate  Eckart  (1951)  solution  is  more  commonly  used  in 
comparisons  with  field  data  due  to  its  simpler  mathematical  formulation.  In  the  field 
data  analysis,  a  difficulty  in  discriminating  between  leaky  and  edge  wave  modes  for 
small  offshore  distances  arises  because  of  the  similarity  of  their  cross-shore  structures. 
Also,  both  waves  have  maximum  amplitude  at  the  shoreline.  Fig.  2.2  shows  such 
similarities  between  higher  mode  edge  waves  and  reflected  waves. 

In  summary,  the  low  frequency  motions  may  be  classified  as  edge  waves,  leaky 
modes,  forced  waves  or  a  combination  as  shown  in  Fig  2.3,  (after  Symonds,  19S2).  A 
complete  set  of  free  wave  solutions  of  discrete  edge  waves  occurs  for  w"  <  gk  while  a 
continuum  of  leaky  mode  exists  for  co  >  gk  .  It  is  noted  that  the  forced  waves  may 
lie  between  the  edge  wave  modes,  but  they  do  not  satisfy  the  dispersion  relationship 
and  are  weaker  than  the  free  modes  since  edge  waves  may  grow  resononantly  as  they 
are  excited. 

2.  Generation  of  Infragravity  waves 

Munk  (1949)  and  Tucker  (1950)  suggested  that  the  long  period  waves  may  be 
caused  by  an  excess  of  mass  carried  forward  by  groups  of  high  swell.  The  swell  waves 
are  assumed  to  be  destroyed  on  the  beach,  and  the  excess  mass  transport  is  reflected 
back  as  a  free  long  wave.  They  tested  this  hypothesis  by  comparing  offshore  pressure 
records,  delayed  by  the  calculated  travel  time  of  the  incident  wave  envelope.  Tucker 
(1950)  observed  a  negative  maximum  correlation  between  incident  wind  wave  groups 
and  the  infragravity  waves  at  a  time  lag  equivalent  to  the  travel  time  of  a  forced  wave 
moving  with  the  wind  wave  group  velocity  and  traveling  back  upon  refiection  as  a  free 
wave. 

Longuet-Higgins  and  Stewart  (1962)  showed  that  the  variations  in  the 
radiation  stress  due  to  incident  wave  grouping  drive  a  forced  wave  such  that  a 
depression  of  the  mean  water  level  (MWL)  will  occur  under  incident  high  waves,  and  a 
corresponding  rise  in  MWL  occurs  under  low  waves.  This  forced  oscillation  of  second 
order  is  associated  with  groups  until  the  incident  waves  break  onshore  and  release  a 
free  wave  offshore,  which  theoretically  supports  the  above  observation. 

Gallagher  (1971)  suggested  edge  wave  excitation  through  a  non-linear 
interaction  whereby  two  incident  waves  can  transfer  energy  into  an  edge  wave  if  the 
difference  frequency  and  longshore  wavenumber  of  the  two  interacting  wind  waves 
satisfy  the  edge  wave  dispersion  relationship,  equation  2.21  .  This  resonance  has  been 
observed  in  the  laboratory  (Bowen  and  Guza,  1978). 
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Figure  2.2     Modes  of  edge  waves  and  reflected  waves,  (p  is  the  velocity  potential. 
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Incident  waves  may  transfer  energy  to  either  leaky  or  trapped  modes  through 
many  possible  interactions.  One  of  these  is  given  by  Guza  and  Davis  (1974)  where 
energy  is  transferred  into  edge  wave  modes  through  a  non-linear  interaction  involving  a 
normally  incident  standing  wave  on  a  sloping  bottom.  Bowen  and  Guza  (1978)  showed 
energy  transfered  into  leaky  modes  (co~  >  gky  )  by  edge  waves. 

Another  mechanism  for  low  frequency  wave  generation  is  suggested  by 
collision  between  the  stronger  component  of  backwash  and  the  shoreward  moving 
bores  that  results  in  roll  waves  which  last  for  few  seconds.  On  analyzing  swash  data, 
Waddell  (1973)  found  that  the  collision  between  uprush  and  backrush  plays  an 
important  role  in  reducing  the  extent  of  swash  runup  waves.  Furthermore,  he 
suggested  that  water  percolation  into  the  beach  serves  as  a  low-pass  filter  affecting  the 
swash  frequency.  Bradshaw  (1980)  suggested  bore-bore  capture  as  a  mechanism  for 
long  wave  generation  in  the  swash  zone.  Mase  and  Iwagaki  (1984)  showed  a 
considerable  shoreward  decrease  in  the  ratio  between  the  number  of  runup  waves  to 
the  incident  waves. 

Symonds  ct  al.  (1982)  developed  a  2-D  model  for  generating  long  waves  which 
is  forced  by  a  time-varying  break-point  caused  by  the  groupness  of  the  incident  wave 
field.  The  generated  free  long  wave  at  the  group  frequency  has  a  standing  wave 
structure  shoreward  of  the  breakpoint  and  an  outgoing  progressive  wave  structure 
seaward  of  the  breakpoint. 

In  the  above  survey,  both  short  and  long  wave  models  are  described  separately 
showing  different  characteristics.  It  is  of  interest  now  to  combine  the  two  waves  and 
study  the  changes  in  the  incident  wave  (short  wave)  field  due  to  the  presence  of  longer 
period  waves. 
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Figure  2.3     Dispersion  relations  for  infragravity  waves. 
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III.  SHORT  WAVE  MODULATION  DUE  TO  LONG  WAVES 

When  short  period  waves  ride  on  much  longer  waves,  they  arc  modulated  in 
amplitude,  wavelength  and  direction.  Their  amplitude  tends  to  increase  at  the  crests  of 
the  long  waves  and  decrease  at  the  troughs.  The  changes  in  wavenumber  and  direction 
depend  mainly  on  kinematical  considerations.  In  this  chapter,  the  changes  in 
wavelength,  direction  and  amplitude  are  investigated.  A  development  lor  the 
corresponding  changes  in  energy  and  momentum  flux  then  follows. 

A.       CHANGES  IN  THE  WAVE  AMPLITUDE 

The  amplitude  modulation  of  a  short  progressive  wave  riding  on  a  much  longer 
wave  was  pointed  out  by  Unna  (1942,  1947).  Unna  (1947)  showed  a  contraction  of  the 
wavelength  and  an  increase  in  the  amplitude  of  short  waves  at  the  long  wave  crests. 

Longuet-Iliggins  and  Stewart  (1960)  used  a  perturbation  analysis  to  examine  the 
non-linear  interaction  between  short  and  long  progressive  waves  over  a  horizontal 
bottom.   They  described  the  surface  elevation 


=      ag   sin  if/g    +   a^  sin  i^£  (3.1) 


where  \|/  =  (kx-  co  t+  0  ),  is  the  phase  function  and  G  the  phase  shift.  The  subscripts 
s  and  £  refer  to  the  short  and  long  wave  respectively.  Using  Stokes'  method  of 
approximation  to  the  second  order,  they  derived  a  general  form  for  the  modulated 
surface  elevation 


n'       =      a    (1   +P)sin    ^      +   ag  Qcos    ^  (3.2) 


where  P  =  a^  k „  sin  \\L  and  Q  =  -a^k$  cos  \|/o  .  Using  the  first  term  in  the  right  hand 
side  of  equation  3.2,  Longuet-Iliggins  and  Stewart  expressed  the  modulated  short  wave 
amplitude  as 
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a'       =      a    {1    +    a^k^J-cothkjh    +    j  tanh  k  ,,h)  sin    ip£}  (3.3) 


and  a  modulation  in  the  short  wave  wavenumber,  k'  ,  as 


k'       =      k    {1    +    a?k^  cothk^h    sirn^j  (3.4) 


In  shallow  water,  equations  3.3  and  3.4  reduce  to 


*'      =      ^i    +   J  -F  sin'^}  (3.5) 


a2, 
k'       =      k    (1    +   -£  sin^}  (3.6) 


These  amplitude  and  wavenumber  modulations  arc  explained  by  the  work  done  by  the 
long  waves  against  the  momentum  of  the  short  waves.  This  work  is  converted  into 
short  wave  energy  and  therefore  produces  a  steepening  of  the  short  wave  at  the  crests 
of  the  long  waves,  i.e.,  the  energy  is  redistributed  along  the  wavelength  of  the  long 
wave.  Longuct-1  Iiggins  and  Stewart  (I960)  came  up  with  equations  similar  to  3.3  and 
3.4  for  a  progressive  short  wave  superposed  on  a  long  period  standing  wave  in  which  a> 
was  described  as  2ui  . 

1.  Determination  of  the  Modulated  Wave  Profile 

It  is  of  interest  now  to  show  that  the  refractive  changes  in  amplitude  and 
wavenumber  due  to  the  non-linear  interaction  may  result  in  a  modulated  short  wave 
expressed  as  an  infinite  sum  of  components  at  different  frequencies.  Then,  the 
properties  of  the  modulated  wave  are  investigated  to  provide  a  better  understanding  of 
the  modulation  process.  The  procedure  is  to  let  the  modulated  short  wave,  designated 
bv  primes,  propagate  in  the  same  wave  direction  and  satisfy 

n'       =      a'    sin  (k'x    -    w't)  (3.7) 
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Substituting  equations  3.3  and  3.4  into  equation  3.7,  and  letting  o    =  (0    ,  yields 


n  '       =      a'    sin  (y      +  mk  sin^£)  (3.8) 


where  nij,    =    a^  k^  ks   x  coth  kp  h    is  defined  as  a   non-dimensional  waVenumber 
(wavenumber  modulation  factor).    Equation  3.8  is  a  general  form  of  amplitude  and 
wavenumber  modulation.    Equation  3.S  reduces  to  the  simple  form  T|    =   a    sin  \\i 
when  m^  <  <   1.    Expanding  equation  3.8  using  trigonometric  identities,  gives 


'       =      a  *  {sin    ty     cos  (m,    sin  >K  )    +   cos    \p     sin  (rru    sin ty . )  }  (3.9) 


Using  Fourier  scries  expansions  to  express  the  sinusoidal  arguments  (Abramovvitz  and 
Stcgum,  1904),  gives 


H'       =      a'         I       J    (mk)sin(i|j      +   nip£)  (3.10) 

n=-°° 

In  the  above  analysis  no  assumptions  have  been  made  regarding  either  applicable 
regions  or  the  initial  properties  of  the  two  interacting  waves.  Equation  3.10  oilers  a 
means  of  expressing  the  modulated  short  waves  as  a  combination  of  spectral 
components  based  on  monochromatic  input.  The  short  wave  is  viewed  now  as  a 
carrier  subjected  to  a  simultaneous  amplitude  and  phase  modulation,  (Panter,  1965). 
Changing  the  limits  of  summation,  equation  3.10  becomes 


a'J0(mk)sin    ^  (311) 


+    a'    .    y      J    (m,){sinU    +n\i>    )    +  (-1)  nsin  (^    -n,j,  Q  )  } 
**  i      n      K.  sjc  s        * 

n=l 
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The  first  term  has  the  form  of  the  solution  for  waves  at  the  incident  wave  frequency 
propagating  on  a  sloping  bottom  but  with  a  different  Bessel  argument  (sec  equation 
2.16  ).  The  remaining  terms  represent  the  side  bands  at  sum  and  difference  frequencies 
(\|/s  ±  n\\h)  where  they  are  phase  locked.  This  solution  is  dependent  on  the  argument 
itik  and  converges  rapidly  as  n  increases. 

In  shallow  water,  the  Bessel  function  argument  m^  reduces  to  2  (a^  /h)  ks  x 
for  the  case  of  progressive  short  waves  superposed  on  long  standing  waves.  The 
argument  mj.  is  found  to  increase  onshore  (Fig  3.1),  where  the  computed  values  of  m^ 
are  evaluated  until  closer  to  the  shore  (x  =  1  cm)  before  mj.  reaches  so  .  Due  to  the 
onshore  increase  of  m^  ,  .ln  (nu  )  decreases  onshore,  oscillating  rapidly  and  reaches 
zero  at  the  shoreline  (Fig  3.2).  Note  the  difference  in  behaviour  of  the  Bessel  functions 
in  equation  3.11  and  2.15,  where  the  Bessel  argument  allows  the  amplitude  of  the  long 
standing  wave  to  reach  a  maximum  at  the  shoreline.  For  breaking  short  waves  on  a 
long  period  standing  wave,  the  depth  limited  breakers  is  utilized  to  describe  a    ( =  y  h) 

in  equation  3.5  .    Substituting  a'  into  3.11,  yields 

3  ai 
n'     =    YhU0(r\)sin  (|/g  +  j  -^  J^m^cos  ipg] 

+  YhiJn(r\){sin(;|;s +rt^)    -    (-l)nsin(i^s  -nji£) } 


(3.12) 


-H(Jn-l(rV  -Jn+1(V){cos(^s+n{,£)  -  (-l)ncos(Vn^)}] 
Two  distinct  features  of  this  solution  arc  demonstrated  in  the  amplitude  squared 
spectra  (Fig  3.3).  First,  symmetric  side  bands  are  generated  about  the  carrier  (Fig  3.3 
a).  The  second  feature  is  the  transfer  of  short  wave  energy  to  the  upper  and  lower  side 
bands  as  the  depth  decreases  (Fig  3.3  b),  which  causes  broadening  of  the  wave  spectra. 
It  is  also  noted  that  the  amplitude  of  the  spectral  components  decrease  shoreward. 
The  instabilities  associated  with  the  side  bands  generation  are  attributed  to  their  spatial 
structure  and  phase  variation.  It  is  shown  that  the  side  bands  energy  is  modified  by 
the  decreasing  depth.  The  modulation  may  contribute  to  the  infragravity  waves  by 
strengthening  the  low  frequency  energy  component  closer  to  the  shoreline,  which  has 
previously  been  observed.  Waddell  (1973)  presents  evidence  of  nonlinear  transfer  of 
energy  from  high  to  low  frequency  due  to  collision,  where  successive  bores  begin  to 
overrun  each  other.  However,  Huntley  and  Bowcn  (1975b)  suggested  that  energy 
transfers  are  due  to  the  interaction  of  many  irregular  spaced  breakers  in  the  surf  zone. 
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Figure  3.1     Wavcnumbcr  modulation  factor  rriu  vs  non-dimensional  distance  x/xh  , 

where  x^  is  the  surf  zone  width.  D 
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Sawaragi  and  Iwata  (1974)  studied  experimentally  the  wave  transformation  after 
breaking  and  observed  energy  transfer  from  the  carrier  into  higher  frequency  waves  by 
an  unknown  mechanism.  The  present  analysis  may  be  used  to  illustrate  some  of  these 
findings,  liquation  3.11  shows  that  the  beach  does  not  eliminate  the  incident  wave 
frequency  and  that  the  modulated  wave  energy  components  and  the  incident  wave 
frequency  vary  spatially. 

B.       CHANGES  IN  THE  WAVE  KINEMATICS 

1.  Background  on  Wave  Kinematics 

At  any  instant  of  time,  a  wave  front  is  defined  by  the  phase  function  \\i  = 
constant,  which  is  the  equation  of  a  family  of  parallel  planes  with  normal  vector  k.  As 
time  increases,  these  planes  move  with  the  phase  speed  in  the  direction  k. 

In  general,  the  dispersion  of  the  wave  motion  may  be  written 


to 


=      w(k,h) 


(3.13) 


where  the  frequency  varies  with  both  the  wavenumbcr  and  the  local  water  depth.  The 
local  wavenumbcr  and  frequency  are  defined  by  the  gradient  of  the  slowly  varying 
phase  function  in  space  and  time: 


k      =    Vij;  (3.14) 


03  = 


3t  (3.15) 


Since  the  curl  of  a  gradient  is  zero,  it  follows  immediately  from  equation  3.14  that  the 
local  wavenumbcr  vector  in  space  is  irrotational,  i.e., 
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V   xk      =      0  (3-16) 


Cross  differentiating  equations  3.14  and  3.15  to  eliminate  V|/  yields 


l£  ♦    V*      =      0  (3-17) 


Equation  3.17  is  known  as  the  kincmatical  conservation  of  wave  density,  where  the  rate 
of  change  of  wavcnunibcr  is  balanced  by  the  convergence  of  the  frequency.  Sometimes 
3.17  is  also  called  the  "conservation  of  crests  equation,"  since  the  crests  arc  neither 
created  nor  destroyed  and  their  total  number  must  be  conserved. 

Rewriting  equation  3.17  in  tensor  notation  and  expanding  the  convergence 
term  using  the  dispersion  relationship  in  equation  3.13  gives 


3ki   +   Jb.3)    +  !tt  -!*L     =      0    .      i-  1.2  O-IS) 


Equation  3.18  can  be  rewritten  as 


oK  .  ok-  _  ^ , 

i.  +    r        -   +    —  -IS-      =0  i    =    1    2 

3t  g.     3x.     +    3h    3xi  U    '  X         L,Z 


(3.19) 


where  d  k-  jd  \{  -  d  kj  jd  X:  from  equation  3.16  ,  and  the  group  velocity  is  defined  by 


3  ui 
q.      =      3kT  (3-20) 
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2.  Kinematics  of  waves  on  current 

In  studying  the  problem  of  short  gravity  waves  superposed  on  waves  of  a 
much  longer  wavelength  and  period,  the  long  wave  may  be  modeled  as  a  slowly  varying 
current  li(x,t),  (e.g.  Longuct-IIiggins  and  Stewart,  1960,  1961;  Garrett  and  Smith, 
1976).  A  coordinate  system  that  moves  with  the  current  velocity  is  chosen,  by  which 
the  observed  frequency  to   ,  passing  a  fixed  point,  is  given  by 


u)'       =      u>     +   k      -    Uz  (3.21) 


where  0>s  is  the  intrinsic  frequency,  (Urscll,  1960;  Whitham,  1960).  The  subscripts  s 
and  1  refer  to  short  and  long  wave  properties.  The  second  term  on  the  right  hand  side 
of  3.21  is  due  to  the  relative  motion  of  short  waves  with  respect  to  the  chosen 
coordinate  and  is  known  as  the  doppler  shift  effect. 

The  short  wave  will  be  modulated  by  the  presence  of  the  current,  and 
therefore  the  modulated  short  wave  (denoted  by  primes)  is  assumed  to  satisfy  equation 
3.17,  i.e., 


^~-  +   Vco'       =      0  (3.22) 


where  k'  ,  the  observed  (modulated)  wavenumbcr,  is  unknown  and  assumed  unsteady 
The  observed  frequency  0)   is  given  by  equation  3.21  .    Ilxpanding  equation  3.22  gives 


!^l  +  c      -3ki   ,   a«u    ah    ,   3kJUJ     =     o  <3-23) 

8t  q  .     3  X  .  dhdX.  d  X  . 

3         3  ii 

where  the  quantities  k-  ,  eg  and  to  arc  for  the  short  waves,  and  U  is  the  current  (long 
waves). 

Assuming  a  steady  intrinsic  short  wave  train;  i.e.,  d  ks  /  3  t  =  0,  equation  3.23 
may  be  simplified  by  the  use  of  equation  3.19  to  give 
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.,  ,  ok      u. 

5k !  s .     1 


L  +  ,1  =      0  0-24) 


3  t  aX. 


which  describes  the  kinematic  conservation  of  the  modulated  wave  density. 

3.  Applications 

The  following  analyses  demonstrate  the  efficiency  and  simplicity  of  equation 
3.24  in  examining  the  wavenumber  changes  of  short  waves  riding  on  a  current  of 
longer  period.  Two  different  bottom  topographies  are  considered:  flat  bottom  and 
plane  sloping  bottom.    Different  current  models  arc  applied. 

a.    Waves  on  a  Current  in  Water  of  Constant  Depth 

The  current  will  be  modeled  as  a  progressive  or  standing  long  period  wave. 

i.  Current  as  a  progressive  long  wave 

The  short  wave  is  assumed  to  propagate  over  a  fiat  bottom  and  to  be 
parallel  to  the  long  wave  (collincar)  in  the  X-dircction.  The  current  U  is  described  by 
the  orbital  velocity  of  a  progressive  long  period  wave, 


cosh   kn  (h+z) 
ui      =      arJi        sinhVh        sin  **  <3^> 


which  corresponds  to  a  free  surface  elevation,  T\~  =    a,  sin  vj/»  .    Substituting  into 
equation  3.24  gives 


3k 


(3.26) 


.     .   cosh  Mhj-z)   ■   ,  -, 

t  +  ks  [a^Ski         sinh  kh    Sin  V 
x  «■ 


Equation  3.26  is  easily  integrated  over  time  to  give 


k'    -    k      k    a      cosh    ¥h  +  z)    cos    i>n       =      Constant  (3.27) 

k  s      11         sinh    kh  *l 

x  X 


43 


where  the  constant  is  evaluated  at  U  =  0,  so  that  k'=  constant  =  ks  ,i.e.  the 
wavenumber  will  not  change  if  there  is  no  current.  Then,  the  general  solution  is 
obtained  such  that 


k'       =      ks     U    +    ^T1  (3.28) 

x  *- 


which  describes  the  wavenumber  changes  due  to  the  non-linear  interaction  between  the 
short  wave  and  a  progressive  long  wave  using  a  kinematic  approach.  Evaluating 
equation  3.25  at  the  free  surface,  z  =  0,  equation  3.28  may  be  written 

k'  =  ks  (  1  +  a£  k£Coth  k^h  sin  \\i»  ) 

It  is  interesting  that  the  above  result,  obtained  from  the  first  principles  of  kinematical 
conservation,  is  the  same  equation  as  (3.4  )  obtained  by  Longuet-lliggins  and  Stewart 
(1900)  using  a  lengthy  second  order  perturbation  analysis  scheme. 

ii.  Current  as  a  long  standing  wave 

For  the  flat  bottom  case,  equation  3.24  reduces  to 


i]Sl   +    k         ^±      =       0  (3.29) 

3t  si    oXi 


The  current  U  is  described  by  the  orbital  velocity  of  a  linear  standing  wave 


cosh  k  (h+z) 
U  -   V<>ke   sinh  k.h Sin  V  Sin  V  (3-30) 


which   varies    only   in   the   x-direction.     Applying   equations   3.29   for   both   x   and  y 
directions  yields 


3k'  •  cosh    k0  (h+z) 

x         ,  I 


-  +    k        a„knu>      ■ — r— -, — r- cos    knx    sin   oj    t      =      0       (3  31) 

s         III         sinh    k,h  I  I  lJ    1J 
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3t 

x  * 


5k' 

"5?     =      ° 

0  t 

Integrating  over  time  gives 


I 


cosh    k    (h+z) 

kv  "  kc  ko   ^^u  v  u "H  =   Constant 

x     s   x,   sinn  k .  n     Is 
x  i  I 


k'   =   k    =   Constant 

y     y 


(3.32) 


Equation  3.32  indicates  that  the  wavenumber  component  alongshore  is  always  constant 

since  no  longshore  component  of  the  current  exists. 

The  constant  in  (3.32)  can  be  evaluated  at  any  node  where  \\»  =   0.    At 

these  locations,  it  is  assumed  that  the  wavenumber  does  not  chance,  i.e.  k'     =   k„    . 

x 
Substituting  into  equation  3.32,  a  general  solution  is  obtained: 


k'       =      k       {1    +   a0k    coth    k„h   cos    knx   cos    w . t }     , 

x  s  I    'l  I  %  i 

x  (3.33) 

|     k'       =      k  =      Constant    . 

V      y  sy 

Therefore  the  wavenumber  component  in  the  x  direction  is  always  modulated  in  time 
and  space,  while  the  longshore  wavenumber  component  is  always  constant.  As  a 
result,  the  wavenumber  vector  changes  in  time  and  space;  i.e.  the  waves  arc  refracted. 

It  is  clear  that  if  the  waves  are  parallel  to  the  current,  i.e.  collinear  in 
x-dircction,  then  k     =  0  and  the  solution  reduces  to  the  first  equation  in  3.33  . 

b.    Waves  on  C uncut  in  Water  of  Variable  Depth 

A  reasonable  model  to  describe  the  current  on  a  sloping  bottom  is  the 
orbital  velocity  of  either  a  cross-shore  long  standing  wave  or  edge  waves.  Collinear 
and  oblique  incident  waves  on  a  current  are  included.  Let  us  assume  a  plane  sloping 
bottom  for  simplicity  of  analysis. 

i.  Current  as  a  cross-shore  long  standing  wave 

On  a  sloping  bottom,  the  horizontal  velocity  of  a  standing  wave  is 
described  by  equation  2.17  .  Recalling  the  equation  for  wave  density  conservation 
(3.24), 
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3ks 
3Ki    +    k         J£   +    u   -4i      =      0  (3.34) 

~oT   +    Ksi    9:<i 


1 


which  requires  specifying  3  V/3  x  and  3  k     (d  x-  .   The  gradient  o[  the  current  velocity 

bi         * 
becomes 


3U  fU  1    3n2-i 

IS     =      "{x   +   h  -Jt]  (3-35) 


Using  the  dispersion  relationship  in  shallow  water  to  differentiate  the  wavenumber 
components  gives 


dk 

_*      =      -   k      (tan    ad^^}  (3-36) 


and  3  ky  /  5  y  =  0,  since  the  dispersion  relationship  is  independent  of  y  as  the  depth 
contours  arc  assumed  straight  and  parallel.  Applying  Snell's  law  to  describe  the 
change  in  the  angle  of  approach  with  x,  3  a  j  3  x  gives 


dx"     =      2h    tan    a    tan    6  (3-37) 


Substituting  into  3.36,  gives 


dk  k  0 

sx  sx  ,         2    . 

IT      =      "    -2x~(seC   a)  (3-38) 


Substituting  into  3.34  using  3.35  and  3.38 

2 

3kx                         'n          1    3n!>                   kqv    sec    a 
Htksl-J-Sirl-»^ 1       -      °  (3.39) 
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3k' 

a  t 


from  which  d  k    jd  y  =  0  and,  therefore,  k'     =  ky  =  constant.    Integrating  equation 
3.39  over  tinic  gives 


k  t  2  n« 

k« 22.  /      u   dt[l    +        ^    a]    -    k        -£     =      Constant  n  dm 

x  x      ;  2  sx      h  (J.4U) 


since 


J  U  dt  =  -(1/(0-  )  d  U/  d  t  +  constant 
Therefore,  equation  3.40  becomes 


.  ^sx  |0U   +  sec!,,    .   k^  ^     .     Constant  (34|) 


The  constant  in  3.41  can  be  evaluated  at  x  =  °o  where  current  vanishes,  i.e.  r\  >  =  0 
and  U  =  0,  which  gives  k'  =  constant  =  k„  .  Then,  for  a  steady  oblique  wave 
train,  the  general  solution  is  given  by 


kx      =      ks    U    +    ~h 2(1    +   — 2~" }     8t]  (3.42) 

x  xw„ 


If  the  waves  and  the  current  are  collinear  (a  =  0),  equation  3.42  reduces  to 

I     kx      =      ks    ll    +   I  "    2    "  It1 

J  x  xoi^  (3.43) 

k'       =      k         =      Constant 

y  y 

This   equation   shows   the   dependence   of  the  wavenumber  changes   on   the   surface 
elevation  and  on  the  acceleration  of  the  current. 
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ii.  Current  as  an  Hdge  Wave 

In  the  3-D  model  the  long  wave  is  described  by  an  edge  wave.  A  zero 
mode  edge  wave  is  used  to  simplify  the  analysis  where  the  Lagucrre  polynomial  equals 
unity.   The  horizontal  velocity  components  are  given  as 


a   qk  -k    X 

U      =      _JZ_Y±  e      y    sin(k      y    -    w.t)  (3.44) 

w„  y0  *- 


-k      X 
V      =      -= — ^  e        l      cos(k,   y   -   aj0t 


oj 


I  2 I 


Yo  l  (3-45) 


A  steady  short  wave  train  riding  on  an  edge  wave  will  satisfy  the  wave  density 
conservation  equation  (3.24  ).  Substituting  into  equation  3.24,  using  equations  3.38 
and  3.39,  gives 


k 
kx    .,„.   „    +  4^  sec2u]   =   0  (3.46) 

x  Ji 


■  -  U[k   k    ■  0w 
3t     l  sv  y. 


3k  ' 

— Y-   +  k   k0   U   =   0 
3t     s  I  (3.47) 


Integrating  both  3.46  and  3.47  with  respect  to  time  gives 


k '    +    k       [( — £)v   +    S^C    -    V]       =      Constant 
x  sx         w£  2XoJ£ 

The  constant  in  3.48,  can  be  evaluated  when  x  -^  oo  where  both  U  and  V  go  to  zero 
and  k'x  =  ks    .    This  leads  to 
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x  I  i 

y  y 


ky      "      ksy{1    +    C£}  (3.50) 


Since  edge  waves  are  progressive  in  the  alongshore  direction,  equation  3.50  is  similar  to 
3.2S  where  the  current  was  modeled  as  a  progressive  wave.  The  changes  in  the 
wavenumber  and  direction  of  short  waves  riding  on  zero  mode  edge  waves  are 
demonstrated  in  Fig  3.4  .  The  modulation  in  the  wavenumber  increases  onshore 
where  the  wavelength  goes  to  zero.  Increasing  the  angle  of  wave  incidence  is  shown 
to  increase  the  wavenumber  modulation. 


C.       CHANGES  IN  WAVE  ENERGY  AND  MOMENTUM  FLUX 

1.  Background 

A  wave  train  composed  of  a  group  of  individual  waves  propagates  with  the 
group  velocity  at  which  the  total  energy  is  transmitted.  Energy  is  usually  expressed  in 
terms  of  average  energy  over  a  complete  wavelength  per  unit  surface  area.  The  total 
energy  of  any  wave  system  is  the  sum  of  both  the  potential  and  kinetic  energy. 

The  average  potential  energy  E    per  unit  surface  area  is  given  by 

t+T         x+L      9 

ed     =     \&l  I         n2  dxdt  (JiS|) 


P 


t  X 


The  average  kinetic  energy  Et.  per  unit  surface  area  is  estimated  by 


t+T         x  +  L  n=0       9       ?       2  -. 

1    _£_       r  r  r         (U    +v   +w    )dz  dxdt  (3.52) 

Ek      "       2    2TL    J  xj  _h; 


For  linear  progressive  surface  gravity  waves,  the  average  total  energy  E-j-  is 
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Figure  3.4     Spatial  changes  in  the  wavenumber  of  short  period  waves 
due  to  the  presence  of  zero  mode  edge  wave. 
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1  2 

ET      =       2    pg    a  ^•53^ 


which  is  equally  partitioned  between  the  potential  and  the  kinetic  energy. 

a.    First  Order  Sea-Swell  Energy  Modulation 

Short  period  waves  riding  on  much  longer  wave  can  be  examined  using  the 
principle  of  superposition.  First  order  modulations  are  examined  using  the  linear 
description  of  both  waves.  The  linear  theory  is  recognized  as  a  reasonable 
approximation  for  its  successful  use  in  describing  many  observed  phenomena  even 
within  the  surf  zone,  (Bowen  et  al.,  1968;  Guza  and  Thornton,  1980).  The  total  surface 
elevation  may  then  be  expressed  as 


nT     =      ns   +    n^  (3.54) 


where  the  subscripts  s  and  1  refer  to  sea-swell  waves  and  long  period  waves 
(infragravity  waves)  respectively.  Inside  the  surf  zone,  short  waves  break,  dissipating 
their  energy  onshore.  The  depth-controlled  wave  breaking  model  is  adopted  (equation 
2.11)  to  describe  first  order  short  wave  dissipation  with  zero  amplitude  at  the 
intersection  of  the  mean  water  level  (MWL)  with  the  beach.  The  depth  in  equation 
2.11  is  modified  to  include  the  slowly  varying  depth  changes  due  to  the  presence  of  the 
long  wave, 


a        =      y  (h  +  T)^)  (3-55) 


where  r| « is  the  surface  elevation  of  low  frequency  wave  motion  on  a  sloping  beach  and 
is  described  by  cither  a  long  standing  wave  or  an  edge  wave.  The  objective  of  this 
approach  is  to  examine  a  first  order  energy  perturbation  (second  order  in  amplitude) 
due  to  the  presence  of  long  waves.  In  the  2-D  model,  the  long  wave  is  modeled  as  a 
long  standing  wave.    The  potential  energy  for  the  combined  short  and  long  waves  is 
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calculated  by  substituting  the  combined  surface  elevation  given  by  equation  3.54  into 
equation  3.51  .  The  kinetic  energy  (equation  3.52)  is  similarly  calculated  using  the 
total  horizontal  velocity, 


uT      =      us    +   u£  (3.56) 


where 


u        =      v  (h    +    r,    )    /g/h   cos   u    t  (3.57) 

s  v,  » 


and  Ui  is  given  by  2.17  .  It  is  assumed  that  the  vertical  velocities  in  shallow  water  arc 
small  compared  with  the  horizontal  velocities,  i.e.,  w  <  <  u  .  The  total  energy  lip 
averaged  over  short  and  long  wave  periods  is  the  sum  of  both  the  potential  and  kinetic 
energy  and  is  given  by 

Em      =      \  pg[(,h)2|  +   \  aj^jjjcx)     +    Ji<x)>l 


T 


«-B 


1      2    2    T2fvi 

4  y  ai  Jo(x) 


(3.58) 


(ws+a)£) 


+    T  Y     ao    Jn(x) I  J 


us"^ 


4    Y     "£    ~0 

The  total  averaged  energy  in  equation  3.58  is  represented  by  four  terms  at  frequency 
bands  tos  ,o)j  ,  ((i)  +  0)i  )  and  (o)s  -  (Oj  ).  The  first  term  in  right  hand  side  is  identified 
as  the  short  wave  contribution  if  the  short  waves  were  considered  alone.  Similarly,  the 
second  term  represents  the  contribution  from  the  long  waves  alone.  The  third  term  is 
due  to  the  nonlinear  interaction  between  the  short  wave  and  the  long  wave  and  is 
represented  by  two  side  bands  at  the  sum  and  difference  between  short  and  long  wave 
frequencies. 
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b.    Model  Parameters  and  Preliminary  Results 

Energy  calculations  for  the  combined  short  and  long  waves  requires 
specifying  the  two  parameters  y  and  aj  in  equation  3.58  .  The  breaker  index  y  ,  which 
relates  the  breaking  wave  amplitude  to  the  local  water  depth,  is  taken  to  be  .20  when 
using  the  root  mean  square  wave  height  (Thornton  and  Guza,  1982).  The  second 
parameter,  the  long  standing  wave  amplitude,  is  defined  as  the  vertical  excursion  of  the 
low  frequency  swash  oscillation  at  the  shoreline.  Earlier  work  by  Munk  (1949)  and 
Tucker  (1950)  showed  the  infragravity  wave  height  to  be  10  percent  of  the  incident 
wave  height,  with  both  quantities  measured  in  15  m  water  depth.  Goda  (1975) 
measured  infragravity  wave  heights  in  1  m  water  depth  and  found  them  to  be  20-40 
percent  of  the  olfshorc  incident  wave  height.  Guza  and  Thornton  (1982)  measured 
infragravity  run-up  height  in  Torrey  Pines  to  be  70  percent  of  the  wind  wave  heights 
measured  in  10  m  water  depth.  Holman  and  Bowen  (1984)  measured  the  significant 
runup  height,  defined  as  the  significant  vertical  excursion  of  water  level  at  the 
shoreline,  and  found  it  to  be  60  percent  of  the  incident  significant  wave  height  at  the 
breakerline.  Guza  and  Thornton  (1985a)  analyzed  data  acquired  from  three  different 
experiments  and  found  that  the  significant  vertical  swash  excursion  on  the  average  was 
equal  to  the  incident  significant  wind  wave  height  measured  in  10  m  depth 


RV      =      H    .  (3.59) 

sig  v        ' 


Therefore,  the  amplitude  of  the  long  wave  a^is  taken  as  .50  II  •  .  This  result  will  be 
used  here  since  equation  3.59  is  weighted  over  three  different  beaches  and  incorporates 
the  Santa  Barbara  data,  which  are  considered  in  the  present  analysis. 

Data  acquired  at  the  Lcadbettcr  Beach  experiment  in  Santa  Barbara  during 
1980  arc  utilized  to  specify  the  parameters  of  the  energy  model.  The  best  fit  to  the 
beach  profile  assuming  a  plane  sloping  beach  is  found  to  be  1:25  with  a  surf  zone  width 
of  50  m.  The  characteristic  periods  for  both  the  short  and  long  waves  as  determined 
from  the  spectral  analysis  of  current  meter  records  arc  14  and  85  seconds  respectively. 
These  arc  typical  values  for  the  data  acquired  on  February  4  (narrow  banded  waves). 
The  model  results  arc  shown  in  Fig  3.5  .  The  potential  energy  for  long  waves  is  shown 
to  reach  a  maximum  at  the  shoreline  with  a  magnitude  comparable  to  the  maximum 
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short  wave  potential  energy  at  the  brcakcrline  (Fig  3.5  a).  The  kinetic  energy  for  both 
short  and  long  waves  decays  shoreward  (Fig  3.5  b).  The  total  averaged  energy  decays 
onshore,  with  short  waves  dissipating  all  their  energy  at  the  shoreline  where  the  long 
waves  have  their  maximum  energy  (Fig3.5  c).  The  two  side  band  components  show  a 
steady  growth  (order  of  15  percent)  in  total  energy  in  the  vicinity  of  the  shoreline. 

2.  Mass  and  Momentum  Fluxes 
a.    Background 

The  changes  in  the  mass  and  the  momentum  fluxes  are  required  to  specify 
the  driving  forces  for  predicting  the  surf  zone  dynamics.  Conservation  equations 
developed  by  Phillips  (1966)  arc  followed,  in  which  the  terms  representing  the  mean 
and  fluctuating  quantities  have  been  separated.  These  equations  are  applicable  to 
wave  motion  as  well  as  general  turbulent  motion.  No  restrictions  arc  placed  on  the 
wave  slopes  and/or  amplitudes  through  the  use  of  these  equations. 

The  conservation  of  total  mass  per  unit  area  can  be  expressed 


2|£  +   -r±     =      0    ,       i   -    1,2  (3.60) 

d  t  dX 


where  i,j  refer  to  horizontal  coordinates,  D  is  the  total  averaged  depth  of  water,  which 
may  include  the  wave  setup,  and  M:  is  the  total  mass  flux  defined  as 


M        =  /       p  u  •  dz    ,       i   -    1,2  (3.61) 

1  -h  X 


where  Uj  is  the  total  horizontal  velocity  that  contains  a  mean  quantity  and  a 
fluctuating  quantity  due  to  the  waves.  The  overbar  denotes  ,  at  this  stage,  averaging  in 
time  over  the  short  wave  periods,  but  not  so  long  as  to  exclude  the  long  waves. 

The  momentum  flux  equations  arc  derived  from  the  horizontal  momentum 
equations  integrated  over  depth,  averaged  over  time  and  using  the  free  surface  and 
bottom  boundary  conditions: 
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Iigure  3.5     Theoretical  wave  model  to  predict  energy  inside  the  surf  zone, 
(a)potcntial  energy  (b)kinctic  energy  (c)to'tal  energy. 
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^   +   MU.M.+S..}       =      "PgD^+R,  (3-62) 

3t         3x.       i    ]         13  ^      3x.  1 


The  first  term  on  the  right  hand  side  is  the  change  in  the  mean  hydrostatic  pressure 
due  to  the  slope  of  the  mean  free  water  surface,  and  R:  is  the  mean  shear  stress  term. 
The  first  term  on  the  left  hand  side  is  the  local  change  in  horizontal  momentum  flux; 
the  second  term  contains  the  mean  momentum  flux  U-  M:  and  the  excess  momentum 
due  to  the  wave  presence  S--  ,  which  was  termed  the  "radiation  stress"  tensor  by 
Longuet-IIiggins  and  Stewart  (1961,  1962).  The  advantage  of  using  the  "radiation 
stress"  technique  to  solve  physical  problems  is  that  the  second  order  effects  are 
obtained  using  first  order  wave  theory,  and  therefore  some  non-linear  wave  properties, 
such  as  setup  and  surf  beat,  can  be  explained  easily.  The  radiation  stress  for  a 
progressive  linear  waves  is  given  by 


Sij      "  /     <-uiuj    +    P6ij)dz    "   2^D  «ij    "    To1  (3.63) 


where  8  is  the  Kroneckcr  delta.  The  third  term  can  be  neglected  since  it  is  of  higher 
order  in  deep  water  and  has  a  negligible  contribution  in  shallow  water.  Evaluating  all 
terms  in  equation  3.63,  the  excess  momentum  tensor  for  a  linear  progressive  wave  is 


F1  c        1^     Jc  d 

S..       =      -?■   -V-   +    \   E[2    -£   -    Ufi.j     ,       i   -    1,2  (3.64) 


where  kx  ,k     arc  the  wavenumber  components  and  are  used  in  the  expression  to 
indicate  the  directions  kY  /k  =  cos  a  and  kv  /k  =  sin  a  . 

These  expressions  arc  also  a  good  first  approximation  for  slowly  varying 
depth  if  a  and  k  are  interpreted  as  the  local  values  corresponding  to  the  local  depth  h 
(Mci,  1983).  In  shallow  water,  the  group  velocity  c  equals  the  phase  speed  c  and  the 
radiation  stress  components  arc  approximated  by 
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=      |[2    cos2a   +11 

2  (3.65) 


=      | [2    sin  a  +11  (3.66) 


(3.67) 
c  E  sina  cos    a 


As  expected,  for  a  standing  wave  on  a  flat  bottom,  the  radiation  stress  components  are 
exactly  twice  the  values  of  a  progressive  wave  (Longuet-lliggins  and  Stewart,  1960). 

b.    Changes  in  The  Radiation  Stresses  Due  to  the  Long  Waves 

The  objective  now  is  to  rederive  the  excess  momentum  due  to  combined 
short  and  long  waves  to  develop  explicit  formulae.  The  short  waves  are  assumed  to 
arrive  obliquely  at  an  angle  a  ,  measured  counter-clockwise  between  the  wave  ray  and 
the  positive  X-axis.  The  long  waves  arc  modeled  as  a  normally  incident  standing  wave 
on  a  plane  sloping  bottom.  The  procedure  of  calculating  the  momentum  flux 
components  due  to  the  presence  of  both  the  short  and  long  waves  is  to  substitute  into 
the  general  formulation  3.63,  which  is  applicable  to  all  kinds  of  steady  and  unsteady 
flow.  In  the  following  development,  the  terms  second  order  in  amplitude  (first  order 
in  energy  and  momentum)  are  retained,  and  all  the  higher  order  terms  are  neglected. 
Then,  the  resulting  expressions  are  averaged  over  the  long  wave  period  to  eliminate  the 
time  dependance  of  the  obtained  formulae. 

The  radiation  stress  component  Sxx  .defined  as  the  flux  of  onshore  directed 
momentum,  is  given  from  equation  3.63  as 


0 !  2 

Sxx      =  /     (puu  +  P)dz    -    j  pgD  (3.68) 

-h 


The  first  term  under  the  integral  is  the  product  of  the  total  horizontal  wave  velocities, 
which  can  be  obtained  by  summing  equations  2.6  and  2.17  .    The  limits  of  integration 
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are  taken  from  the  bottom  to  the  still  water  level  (SWL)  at  \\  =  0,  which  is  consistent 
with  the  analysis  being  second  order  in  wave  slope.  Integrating  over  depth  and 
averaging  in  time  over  long  wave  periods  yields 


2-r2fv1 
2U2        2,7   JCTa;         T2(vn  )  (3-69) 

dz      =      i-pg(v    h     +a?( 2 +  J1«))) 


/       puuTdz      =      2pglY    n  l'         2 

•h 


The  pressure  term  in  equation  3.68  is  obtained  by  integrating  the  vertical  momentum 
equation  (neglecting  the  vertical  stresses)  from  any  depth  z  to  the  free  surface  (see 
Phillips,  1966),  which  gives 


n  g  n  g  n  2 

P|         «      p{    /      gdz    +  tt        /      wdz    +   ^-      /      uwdz    -    pwr  }    (3J0) 
z  z  z  j    z 


The  first  term  is  the  hydrostatic  pressure.     Integrating  from  the  bottom  to  the  free 
surface  and  averaging  over  short  and  long  wave  periods  gives 


7     7  2 

term    1      =      Ipg   :l~y~  +    (n+h)2    +   \{\   +    l)a2JQ(X)}  (3  ?1) 


The  second  term  is  similarly  evaluated,  but  it  was  found  to  vanish  when  averaged  over 
the  long  wave  period.   The  third  and  fourth  terms  give 


12    2       2  2 

term    3      =  -jpha^  u^  [  Jq  ( x)    -J^(X)]  (3.72) 

+   jpu   tan  3    /gh   a,  JQ  (>;)  J^  (>;) 


(3.73) 


7    7  2 

7      2     2       2            £    0    A    ,  2       HO          i 

term    4      =      ph^^y     (h      +    2 )    +    u£    2 ; 
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The  time  averaged  radiation  stress  S„x  is  obtained  by  adding  equations 
3.69,  3.71,  3.72  and  3.73  : 


S 
xx 


=      ^Pg{3Y2h2    +    (1    +   2£_)aj[j*(x)    +    2a\j\{,)} 
+  ^-pu^a2(tan   3    /gh   JQ(x)J1(X))    +    gP^    h 


(3.74) 


+ 


2 

|a2,h{'-^(jQ(x)+J2(x))    +^2J2(X)} 


The  first  term  on  the  right  hand  side  is  identified  as  the  short  wave  contribution  if 
considered  alone,  which  is  in  agreement  with  equation  3.64  obtained  by 
Longuet-IIiggins  and  Stewart  (1962).  The  second  and  the  third  terms  are  the 
contributions  by  the  long  waves  and  the  side  bands. 

The  transverse  component  of  radiation  stress,  S  Y  is  defined  as  the  onshore 
flux  of  longshore  momentum  across  a  vertical  plane  any  distance  offshore.  This 
component  can  be  calculated  from  the  direct  evaluation  of  the  momentum  flux  as  the 
product  of  the  total  horizontal  velocities,  u-r  and  v-y  , 

ni    

Syx      =      J         UTVTdz  (3-75) 


A    cross-shore    long   wave    is   considered    for   simplicity,    i.e.    no    longshore   velocity 
component,  while  an  oblique  short  wave  is  considered.   Then,  equation  3.75  becomes 


;  p         f       (u     cos  a  +u„)u^  sin  a   dz  (3  76) 

yx  u  s  x.        S  v^-v/ 


-h 


The  transverse  component  S      is  the  driving  force  for  the  longshore  current  inside  the 
surf  zone.    When  considering  the  cross-shore  long  waves,  the  primary  change  in  S 
arises  from  the  ainplitudc  modulation  of  the  short  waves,  which  is  a  function  of  the 
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long  waves.    Integrating  equation  3.76  over  depth  and  averaging  over  the  short  wave 
period  gives  the  time  dependent  form, 

Syx      =      |pg{Y2(h+n,)2)sin   a   cos   a  (j  ??) 


Averaging  over  the  long  wave  period  to  eliminate  the  time  dependency,  the  transverse 
component  Syx   inside  the  surf  zone  is 


c  L       ,    2,  2         12    2.2.    n    ,  (3.78) 

vx      =      2pg    {  2Y    a«Jo^^sin    a    cos    a 


The  first  term  on  the  right  hand  side  is  in  agreement  with  equation  3.67,  while  the 

second  term  can  be  identified  as  side  bands.    In  equation  3.7S  ,  it  is  noted  that  the 

cross-shore  long  wave  by  itself  does  not  contribute  to  S       ,  but  their  interaction  with 

the  short  wave  results  in  the  side  bands  that  increase  the  Svv  term. 

> A 


60 


IV.  DATA  ANALYSIS 

The  analysis  of  field  data  is  performed  to  test  the  hypotheses  that  the  nearshore 
dynamics  are  due  to  both  the  short  and  long  waves.  Outside  the  surf  zone,  waves 
arriving  in  groups  cause  a  change  in  MWL  (Fig.  4.1  )  and  drive  a  forced  long  wave 
(Longuet-Higgins  and  Stewart,  1962).  At  the  breaker  zone,  the  short  waves  break  and 
their  heights  become  depth  dependent,  obeying  the  saturation  curve.  The  water  depth 
is  modified  by  the  presence  of  the  long  waves.  Therefore,  the  short  waves  are  expected 
to  be  modulated  by  the  infragravity  waves  inside  the  surf  zone.  The  experiments  are 
briefly  discussed,  followed  by  analyzed  results. 

A.       DESCRIPTION  OF  THE  FIELD  SITES 

Two  experiments  were  performed  as  a  part  of  the  Nearshore  Sediment  Transport 
Study  (NSTS).  The  objective  of  NSTS  was  to  develop  an  improved  engineering 
formula  to  predict  sediment  transport  on  beaches.  Field  sites  were  selected  for  their 
simple  topography  of  straight  and  parallel  depth  contours,  which  match  the 
assumptions  often  used  for  nearshore  dynamics.  The  first  experiment  was  located  at 
Torrey  Pines  Beach  (T.P.),  San  Diego,  California  and  conducted  for  one  month  in 
November,  1978.  Torrey  Pines  is  a  gently  sloping  beach  (tan  p  =  0.002)  with  nearly 
straight  and  parallel  contours.  The  bottom  is  composed  of  moderately  sorted, 
fine-grain  sand  (mean  diameter  0.15-0.2  mm).  The  waves  break  as  spilling  or  mixed 
spilling-plunging.  A  wide  variety  of  wave  and  weather  conditions  were  encountered 
during  the  experiment,  from  small  to  large  (2  m)  waves,  from  very  narrow-band  swell 
to  wide-band  sea,  from  calm  to  windy  days  (>  10  m/s).  A  total  of  42  sensors  were 
deployed  over  an  area  520  m  long  parallel  to  the  shore  and  up  to  500  m  offshore. 
Surface  elevation  and  horizontal,  orthogonal  velocity  components  were  measured  along 
a  shore-normal  transect  from  offshore  at  the  10  m  depth  contour  to  the  shore. 

Leadbcttcr  Beach,  Santa  Barbara  (S.B.),  California  was  the  second  NSTS 
experimental  site  (February,  1980).  Leadbcttcr  Beach  has  relatively  straight  and 
parallel  nearshore  depth  contours  with  well-sorted  fine  to  medium  size  sand.  The  mean 
nearshore  slope  varied  between  0.017  and  0.05  during  the  experiment,  depending  on  the 
wave  climate.    Well-developed  cusps  occurred  during  the  beginning  of  the  experiment. 
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Figure  4.1     Schematic  diagram  of  short  and  long  wave  models  outside 

and"  inside  the  surf  zone. 
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No  major  offshore  bar  system  was  apparent.  The  shoreline  has  an  unusual  east-west 
orientation  along  a  predominantly  north-south  coast.  The  open  ocean  waves  are 
limited  to  a  narrow  window  of  approach  (±  9°  centered  on  249°)  because  of 
protection  from  Point  Conception  to  the  north  and  the  Channel  Islands  to  the  south. 
The  generally  highly  filtered  ocean  swell  type  waves  from  almost  due  west  must  make  a 
right  angle  turn  to  approach  the  beach  normally.  The  result  is  often  a  narrow  band  (in 
frequency  and  direction)  spectrum  of  waves  approaching  at  large  oblique  angles  to  the 
bottom  contours  in  the  surf  zone  and  driving  a  strong  longshore  current.  Rod  and 
level  surveys  of  the  beach  profile  were  obtained  daily.  An  array  of  current  meters  and 
pressure  sensors  were  deployed  extensively  in  a  similar  manner  as  at  T.P.  Severe 
storms  occurred  during  the  experiment,  causing  significant  changes  in  the  shoreline 
configuration  and  beach  profile.  Despite  the  destruction  of  most  of  the  instruments 
during  the  storm,  a  very  substantial  and  valuable  data  was  recorded  (Gable,  1980). 

Data  were  acquired  at  a  sample  rate  of  64  Hz  for  several  hours  each  day,  then 
low-pass  filtered  and  reduced  to  2  Hz.  The  sensors,  data  acquisition  system  and 
experimental  sites  are  fully  described  by  Gable  (1979,  19S0). 

B.       DATA  ANALYSIS 

Data  acquired  from  T.P.  and  S.B.  have  been  extensively  analyzed  by  different 
investigators.  Incident  wave  analyses  include  shoaling,  setup,  swash  oscillation, 
runup,  velocity  moments,  longshore  current  and  sediment  transport.  Infragravity 
waves  were  studied  for  runup,  identifying  both  the  cross-shore  and  longshore  structures 
and  type  of  the  dominant  motion,  whether  they  are  leaky  and/or  edge  waves  (Guza  and 
Thornton,  1982,1 985a, b;  Oltman-Shay  and  Guza,  1986).  As  mentioned  before,  the 
dynamics  of  the  infragravity  waves  were  not  previously  studied.  Furthermore,  their 
effect  on  the  incident  waves,  especially  inside  the  surf  zone,  is  of  interest  to  understand 
the  nearshore  processes.  It  is  hypothesized  in  this  study  that  the  interaction  between 
the  short  and  infragravity  waves  generates  other  waves  (side  bands)  at  the  sum  and 
difference  frequencies.  These  waves  occur  at  the  neigbourhood  of  the  short  wave 
frequency,  and  they  arc  of  importance  in  driving  longshore  current  and  sediment 
transport  at  the  shoreline.  The  existence  of  both  the  infragravity  waves  and  the  side 
bands  are  investigated  utilizing  records  of  wave  staffs,  current  meters  and  pressure 
sensors.    Coupling  between  the  incident  waves  and  the  infragravity  waves  arc  tested  via 
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spectral  and  cross-correlation  techniques,  and  different  theories  arc  discussed.  Since 
the  present  developed  model  is  monochromatic,  only  days  of  narrow  banded  waves  are 
considered.  Three  days  are  selected  from  each  experiment  (total  of  6  days):  3,  4,  and  6 
February  1980  at  S.B.  and  10,  20,  and  21  November  1978  at  T.P. 

1.  Spectral  Analysis 

Velocity  and  pressure  time  scries  of  68-minute  length  are  used  to  compute 
spectra  with  32  degrees  of  freedom  and  a  resolution  of  0.004  Hz.  The  power  spectra 
arc  computed  for  all  the  available  working  sensors  at  each  selected  day.  The  spectra 
are  plotted  on  linear  scale  to  demonstrate  the  relative  importance  of  each  peak  and 
normalized  by  the  variance  (area  under  spectra).  In  Fig  4.2  ,  a  well  defined  swell  peak 
is  seen  at  a  frequency  around  f  =  .06  Hz  in  most  of  the  instrument  locations  during 
10  Nov.,  T.P.,  and  in  Figures  (4.3,  4.4)  at  frequency  f  ~  .07  Hz  in  all  the  analyzed 
instruments  in  3  and  4  Feb.,  S.B.  A  consistent  valley  near  0.05  Ilz  is  present  in  almost 
all  spectra,  separating  the  incident  wave  band  from  the  low  frequency  band.  Within 
the  infragravity  band,  a  pronounced  low  frequency  peak  is  shown  in  most  of  the 
spectra  at  L  =  .017  Hz  for  4  Fcb.S.B.   and  .0195  Ilz  in  10  Nov./f.P. 

'I  he  infragravity  energy  has  been  studied  extensively  since  the  first 
observations  by  Munk  (1949)  and  'fucker  (1950).  Huntley  ct  al.  (1981)  identified  a 
progressive  low  mode  edge  wave  in  the  longshore  direction  that  satisfies  the  dispersion 
relationship  of  edge  waves.  Oltman-Shay  and  Guza  (1986)  extended  the  work  of 
Huntley  et  al.  (1981)  and  were  able  to  determine  the  energy  contained  in  each  edge 
wave  mode.  The  low  mode  edge  waves  (n  ^  2)  were  found  to  dominate  the  longshore 
current  energy  while  leaky  mode  or  high  mode  edge  waves  dominated  the  cross-shore 
direction.  They  attributed  the  daily  changes  in  infragravity  energy  levels  to  the 
variation  in  incident  wave  conditions  and  the  resonant  forcing  of  edge  waves  described 
by  Gallagher  (1971).  Guza  and  Thornton  (1985a)  observed  surf  beat  in  3  field  sites 
including  S.B.  and  T.P.  They  found  a  phase  shift  of  n  /2  between  the  surface  elevation 
and  the  horizontal  velocity  within  low  frequency  motions  (f  ^  .05  Ilz)  which  indicates 
the  predominance  of  the  standing  wave  structure  of  surf  beat.  This  result  agrees  with 
Suhayda  (1974).  Guza  and  Thornton  (1985a)  rejected  the  hypothesis  that  surf  beat  is 
primarily  generated  by  bore-bore  capture  within  the  surf  zone  (as  hypothesized  by 
Bradshaw,  1980  and  others)  because  bore-bore  capture  models  do  not  predict  incoming 
surf  beat  energy  outside  the  surf  zone  and  do  not  describe  the  gradual  increase  in  the 
infragravity  band  onshore. 
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Figure  4.2     Normalized  power  spectra  at  different  locations  in  10  Nov.,  T  .P., 
showing  an  onshore  relative  increase  in  the  inlragravity  energy  level. 
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The  transformation  of  the  incident  wave  spectra  across  the  surf  zone  is  shown 
in  Figures  4.2  and  4.3  .  A  dramatic  shoreward  increase  in  the  low  frequency  band 
energy  is  accompanied  with  rapid  attenuation  in  the  progressive  incident  wave  energy 
due  to  breaking  as  the  depth  decreases.  The  infragravity  energy  can  exceed  the  wind 
wave  energy  by  a  factor  of  2.5  in  the  inner  surf  zone  (Wright  et  al.,  1982).  In  storm 
conditions,  Ilolman  et  al.  (197S)  observed  a  considerable  increase  in  the  infragravity 
energy  levels  responding  to  the  changes  in  incident  wave  conditions.  It  is  clear  from 
Figures  (4.2,  4.3  and  4.4)  that  the  spectra  closer  to  the  shoreline  is  dominated  by  the 
low  frequency  band. 

The  theory  of  wave-wave  interaction  presented  in  Chapter  3,  by  which  side 
band  waves  are  generated  as  a  result  of  the  non-linear  interaction,  can  be  tested  using 
the  available  data.  Flgar  and  Guza  (1985),  on  studying  the  nonlinear  dynamics  of  the 
shoaling  waves  using  bispectral  analysis,  showed  that  the  low  frequency  motion  at  the 
beach  face  is  significantly  non-lincarly  coupled  to  higher  frequency  modes  seaward  the 
surf  zone  around  the  peak  of  the  power  spectra.  They  did  not  resolve  possible  side 
band  generation  in  their  analysis,  choosing  instead  to  obtain  higher  stability  (low 
resolution).  Using  higher  resolution  here,  peaks  at  frequencies  (f  ±  f,  )  are  found  in 
Fig.  4.5  which  are  identified  as  upper  and  lower  side  bands.  This  supports  the 
hypothesis  of  the  short  wave  modulation  derived  in  equation  3.11  .  Evidence  of  side 
bands  is  present  for  almost  all  the  measurements  on  all  days  analyzed,  but  usually  only 
the  upper  side  band  at  frequency  (f  +  Q  )  is  present  (sec  Figures  4.2,  4.3  and  4.4). 
The  side  bands  arc  not  always  statistically  significant,  however,  since  at  least  the  upper 
side  band  is  present  for  all  narrow  band  days.  Hence,  it  is  concluded  that  the  side 
bands  are  real.  It  is  not  clear  why  the  lower  side  band  is  absent.  The  suppression  of 
the  side  band  may  be  strong  lower  frequency  breaking  or  it  may  be  caused  by  the 
seaward  radiation  of  energy  by  the  outgoing  component  of  the  standing  long  waves. 
Flgar  (1985)  showed  that  the  low  frequency  reflected  waves  will  have  phase 
relationships  relative  to  their  incoming  parts  that  tend  to  decrease  bispectral  level. 
Guza  et  al.  (1984)  show  an  increase  in  coherence  between  the  incident  wave  envelope 
and  the  long  standing  waves  when  removing  the  reflected  wave. 

The  time  history  of  energy  transferred  between  swell  peak  and  the  upper  side 
band  is  demonstrated  in  Fig.  4.6  using  almost  a  four-hour  time  scries.  The  long  record 
is  divided  into  64  segments  each  of  17.1  minutes  with  a  50%  overlap  to  calculate  the 
spectra.     Energy  appears  to   be  transferred   from  the  incident  wave  peak  to  higher 
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Fieure  4.3     Normalized  power  spectra  at  different  locations 

in  3  Feb.,  S  .B. 
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Figure  4.4    Normalized  power  spectra  at  different  locations 

in  4  Feb.,  S.1J. 
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Figure  4.5     Normalized  spectra  showing  significant  peaks  at 
frequencies  of  dominant  long  wave  fp,  short  wave  i    and  side  bands. 
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frequency  components.    This  leads  to  a  faster  growth  of  the  upper  side  band  than  for 

the  lower  side  band.    The  incident  wave  peak  drops  to  a  minimum,  strengthening  the 

upper  side  band.    A  similar  mechanism  could  happen  if  lower  side  bands  exist,  which 

may  transfer  energy  to  low  frequency  components. 

Using  the  same  data  set,  Guza  and  Thornton  (19S5b)  computed  the  velocity 

variance  in  20  Nov.,  T.  P.  at  each  instrument  location  across  the  surf  zone  Fig.  4.7  . 

Their  result  show  no  substantial  differences  between  the  observed  and  the  predicted 

values  when  summing  the  potential  and  kinetic  energy  for  short  waves  riding  on  long 

standing  waves  (Fig.  3.5).   The  decrease  in  short  wave  energy  is  almost  balanced  by  the 

increase  in  low  frequency  energy  resulting  in  an  almost  constant  cross-shore  velocity 

variance  (Guza  and  Thornton,  1985b).   Wright  et  al.  (1982)  obtained  a  similar  variance 

distribution  on  Australian  beach.   Although  the  spectral  shape  of  3  Feb  is  similar  to  4 

Feb,  waves  in  the  4   Feb  are  more  energetic  (H  ■      =    80  cm).    Therefore,   strong 

sit, 

non-linearities  arc  expected. 

2.  Cross-Correlation  Analysis 

The  objective  of  this  analysis  is  to  cross-correlate  the  incident  wave  envelope 
with  the  long  waves  to  test  the  hypotheses  that  long  waves  are  forced  outside  the  surf 
zone  and  short  waves  modulated  inside  the  surf  zone.  The  cross-correlation  function, 
Cn  ,  is  a  measure  of  the  degree  of  linear  relationship  (association)  between  two  data 
sets  in  the  time  domain.  Tucker  (1950)  correlated  the  incoming  wave  envelope  and  the 
corresponding  long  wave  observed  by  a  bottom  mounted  pressure  sensor  at  a  distance 
800  m  offshore.  Me  obtained  a  negative  maximum  correlation  corresponding  to  a  5 
minute  time  lag.  The  time  lag  associated  with  the  maximum  correlation  was 
approximately  the  time  required  for  the  incident  wave  group  to  reach  the  shore  with 


velocity  c  and  for  the  associated  long  wave  to  travel  back  with  phase  velocity  V  gh. 
The  negative  correlation  was  described  by  Longuet-lliggins  and  Stewart  (1962)  as  a 
low  frequency  forced  wave  response  due  to  the  incident  wave  groupness  that  has  a 
phase  shift  of  K  with  the  envelope  of  the  high  frequency  waves.  Hence,  under  high 
waves,  a  long  wave  trough  is  expected  and  vicc-vcrsa. 

Iluntly  and  Kim  (1984)  obtained  a  positive  correlation  between  long  waves 
and  the  short  wave  envelope  at  zero  time  lag  from  current  meter  measurements  taken 
outside  the  breaker  line.  Their  results  agree  with  the  forced  wave  model  of 
Longuet-lliggins  and  Stewart  (1962).  For  velocity  measurements,  a  positive  correlation 
is  expected  with  the  x-dircction  positive  offshore. 
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rigure  4.6    Spectral  time  series  of  current  meter  CODx  with  50%  overlapping 
calculated  from  four  hour  record,  (4  Feb,  1980,  Santa  Barbara). 
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Guza,  ct  al.  (1984)  obtained  a  negative  correlation  outside  the  surf  zone  and  a 
positive  correlation  inside  the  surf  zone  when  comparing  the  surface  elevation  of  the 
short  wave  envelope  with  long  waves.  No  explanation  was  given  for  change  in  signs  of 
the  correlation  between  outside  and  inside  the  surf  zone. 

Surface  elevation  time  series  are  inferred  from  current  meter  data  from 
experiments  at  T.l\  and  S.B.  using  the  complex  Fourier  spectra  of  the  horizontal 
velocity  component,  U(f).  Current  meter  records  were  first  linearly  detrended  to 
exclude  the  effects  of  the  rising  and  falling  of  the  tides,  then  high-pass  filtered  using  a 
cut-olf  frequency  of  0.05  1  Iz  to  exclude  surf  beat,  followed  by  a  low-pass  filter  with  a 
high  frequency  cutoff  (0.31  Iz)  to  exclude  capillary  waves.  The  filtering  was 
accomplished  by  Fourier  transforming  the  signals  and  zeroing  the  Fourier  amplitude 
coeilicicnts  in  the  filtered  frequencies.  The  complex  surface  elevation  spectrum,  X(f), 
is  calculated  applying  the  linear  wave  theory  transfer  function,  Ilj-(f) 

x(f)        =       H     (f)U(f) 

f  (4.1) 

Then,  the  complex  surface  elevation  spectrum  is  inverse  transformed  to  obtain  the 
surface  elevation  time  series.  The  entire  68-minute-record  was  transformed  at  one  time 
to  minimize  the  end  effects,  which  result  in  spectral  leakage,  and  to  obtain  maximum 
resolution  for  very  sharp  roll-oir  at  the  filter  cut-offs.  Surface  elevations  were  also 
computed  from  pressure  signals  in  a  similar  manner  by  transforming  the  pressure 
records  using  linear  theory. 

Using  the  surface  elevation  time  scries  obtained  in  this  manner,  the  correlation 
between  the  incident  wave  envelope  and  the  corresponding  low  frequency  motion  is 
calculated.  The  low  frequency  motion  is  obtained  by  low-pass  filtering  the  surface 
elevation  time  series.  The  high  frequency  wave  signal  is  the  residual  after  subtracting 
off  the  low  frequency  energy.  The  short  wave  envelope  is  obtained  by  demeaning  and 
squaring  the  high  frequency  signal,  followed  by  the  same  low-pass  filter  used  before. 
Segments  of  the  time  series  for  the  wave  envelope  and  the  associated  low  frequency 
motion  suggest  positive  correlation  for  the  wave  sensor  inside  the  surf  zone  (Fig.  4.8) 
and  negative  correlation  for  the  current  meter  derived  signal  outside  the  surf  zone  (Fig. 
4.9). 

The  cross  correlation  function  is  calculated  using  the  following  formula 


72 


Figure  4.7     Spatial  distribution  of  band  passed  velocity  variance  in  20  Nov., 
b  Low  is  (0<  1  <0.(b  Hz)  and  high  is  (0  05 <f<  .5  Hz) 

(alter  Guza  and  Ihornton,  1985b). 
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Figurc  4.8     Segment  of  the  time  series  for  the  wave  envelope  (solid  curve)  and  the 
associated  low  frequency  (dotted  curve)  for  wave  staff  W3  8  (inside  the  surf  zone) 

in  21  Nov.,  1978,  Torrey  Pines. 
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in  4  Feb.,  1980,  Santa  Barbara. 
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k    =    0,1,2,  ..  . 

(4.2) 


where  T  =  k  t  is  the  time  lag,  t  is  the  sampling  time  interval,  n  is 
number  of  data  points,  <7j  <T  are  the  standard  deviations  for  the  low  frequency  signal 
and  the  envelope,  and  x,  y  are  the  means  for  each  of  the  time  scries.  The  time  lag  is 
positive  when  the  long  waves  are  in  advance  of  the  wave  envelope,  and  it  is  negative 
when  long  waves  lag  the  wave  envelope. 

A     test     of     significance     (null     hypothesis)     based     on     the     "Fishcr-Z 


Transformation"  is  performed  with  the  statistic  z  =  J  n-3  /2  In  (1+C^  /  1-  Cn  ) 
(Miller  and  Freund,  1977).  At  level  of  significance  a  =  0.05,  the  calculated  z  always 
exceeds  z„/2  (  =  1.96),  since  the  used  record  is  long  (S192  data  points).  Thus,  the  null 
hypothesis  is  rejected  and  the  correlation  is  significantly  different  than  zero. 

The  correlation  function  for  the  wave  stall  W38  (T.P.)  is  shown  in  Fig.  4.10  . 
The  maximum  correlation  is  +0.49,  corresponding  to  a  time  lag  of  -4.0  sec.  The 
function  shows  a  strong  periodicity  at  about  40  sec.  The  correlation  coefficient  is 
calculated  with  95%  confidence  to  be  significantly  different  from  zero.  Fig.  4.1 1  shows 
a  negative  correlation  for  C07x  (S.B.),  indicating  the  characteristics  of  a  forced  wave 
response  outside  the  surf  zone. 

The  maximum  correlation  coefficients  and  the  corresponding  time  lags  at  each 
instrument  location  in  10  and  21  November,  T.P.  arc  given  in  Table  1  and  II  . 
Similarly,  tables  111  and  IV  are  given  for  4  and  6  Fcb.,S.B.  Lower  values  of  Cj^  imply 
that  the  two  random  variables  may  be  independent,  that  there  is  no  linear  relationship 
between  the  two  signals,  suggesting  a  free  long  wave  (Huntly  and  Kim,  1984)  or  the 
possibility  of  some  non-linear  relationship.  The  low  values  of  the  negative  C»  can  be 
attributed  to  the  fact  that  the  n  phase  shift  is  only  a  deep  water  case,  while  the  deepest 
instrument  used  is  3  m.  The  other  reason  could  be  the  significance  of  the  phase-locked 
reflected  low  frequency  waves  decreasing  the  coherence  even  for  low  values  of 
reflection  coefficient  (Flgar,  1985).  Flgar  and  Guza  (1985)  indicated  that  the 
infragravity  modes  do  not  appear  to  be  bound  with  a  fixed  phase  relationship  to  the 
high  frequency  wave  groups,  since  their  biphase  evolves  towards  lower  values  as  waves 
shoal. 

76 


/ 

0 

f 

0 
-8 

0 
-8 

( 

/ 

0 
-9 

oo 

0.0            20.0 

lag   (sec) 

CO 

3: 

1 

E 
0  +-> 

0 

-s 

1 
0 

m4 

0 
•? 

0 

8 

OO'I 

SZO          OS'O          SZ'O          OO'O          STO-        OS'O- 

•^aco  uoi.;e[ajjoD-ssoj3 

sro-    oo-i 

Figure  4.10     The  cross-correlation  function  for  W38,  21  Nov.,  1978, 
Torrey  Pines.   A  positive  maximum  correlation  (.49)  corresponds  to  -4.0  time  lag. 
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Figure  4.1 1     The  cross-correlation  function  for  C07x,  4  Feb.,  1980, 
Santa  Barbara.   A  negative  maximum  correlation  (.29)  corresponds  to  +  3.0  time  lag. 
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It  is  of  interest  to  observe  a  consistent  pattern  of  CR  ,  being  negative  outside 
the  surf  zone  and  positive  closer  to  the  shoreline  (Figures  4.12  and  4.13).  The 
statistically  significant  values  of  CR  in  the  ncarshore  zone  at  zero  time  lag  indicated 
significant  coupling  between  the  high  and  low  frequency  motions.  A  zero  correlation 
in  Fig.  4.12  is  found  to  coincide  with  the  approximate  location  of  the  breaker  line.  The 
zero  correlation  could  be  due  to  the  two  signals  being  non-linearity  related, 
randomness  (noise)  introduced  into  the  wave  signal  by  the  breaking  wave,  a  possibility 
that  the  breakpoint  could  be  at  a  nodal  point  for  the  long  standing  wave,  independence 
of  the  two  signals,  or  a  transition  location  between  forced  and  free  long  waves.  Fig. 
4.13  shows  a  consistent  negative  correlation  outside  the  surf  zone  and  a  positive 
correlation  shoreward  of  the  midsurf  zone.  This  could  mean  that  the  wave  breaking 
docs  not  destroy  completely  the  wave  groupness  until  some  location  inside  the  surf 
zone.  This  location  might  be  a  transition  between  the  forced  and  free  long  waves. 
Flgar  (1985)  showed  that  surf  beat  biphascs  decreases  as  waves  shoal.  On  some  days, 
the  correlation  in  Fig.  4.13  shows  a  rapid  decrease  when  waves  get  closer  to  the 
shoreline.  This  might  be  attributed  to  the  bore-bore  capture  mechanism.  Mase  and 
Iwagaki  (1984)  illustrated  that  a  wave  (bore)  could  be  captured  by  a  successive  bore 
before  the  bore  reaches  a  maximum  runup  height.  Therefore,  the  number  of  runup 
waves  and  frequency  of  runup  will  be  reduced  compared  to  that  of  incident  waves.  The 
ratio  of  number  of  waves  is  given  using  the  surf  similarity  parameter  and  shows  a 
decrease  in  number  of  waves  reaching  a  gentle  slope  beach. 

The  positive  correlations  obtained  inside  the  surf  zone  are  demonstrated  by 
the  simple  model  developed  in  chapter  3  by  linear  superpose  the  two  waves.  The 
surface  elevation  time  series  inside  the  surf  zone  is  simulated  using 

nT      =      Y  (h  +nl)cos    ips    +    r)z  (4.3) 

The  wave  envelope  and  low  frequency  motion  arc  obtained  as  described  above.  A 
positive  correlation  is  always  shown  at  zero  time  lag  (Fig.  4.14  ),  which  agrees  with  the 
observations. 

In  the  following  chapter,  the  theoretical  results  obtained  in  Chapter  3  arc 
applied  to  derive  wave  setup,  longshore  currents  (steady  and  unsteady)  and  sediment 
transport  due  to  the  combined  short  and  long  period  waves. 
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Figure  4.12     Spatial  distribution  of  cross-correlation  coefficients 

snowing  a  dominant  negative  correlation  outside  the  surf /one 

and  a  positive  correlation  inside  the  suil  /.one.    H.  I.,  is  the  1  icakerline. 
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Figure  4.14    Theoretical  model  to  examine  the  positive  correlation 

between  the  wave  envelope  and  the  associated  low  frequcnev  motion  within  the 

surf  zone,  (a)  segment  of  time  scries  (b)  cross-correlation  function. 
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V.  APPLICATIONS 

A.       WAVE  SETUP  AND  SETDOWN 
1.  Background 

In  the  nearshore  zone,  it  is  observed  that  waves  during  shoaling  and  breaking 
produce  a  variation  in  the  mean  water  level  (MWL).  This  variation  in  MWL  is 
considered  a  primary  cause  of  nearshore  circulation,  such  as  rip  currents.  The  change 
in  MWL,  r\  ,  is  characterized  by  two  regions  (a)  a  gradual  depression  of  the  sea  water 
level  starting  from  offshore,  and  reaching  a  maximum  at  the  breakerlinc,  and  (b)  an 
increase  in  the  MWL  shoreward  of  the  breaker  line.  The  depression  is  termed 
"sctdown"  and  the  upward  slope  is  called  "setup." 

To  study  the  wave  setdown  and  setup,  the  X-componcnt  of  the  horizontal 
momentum  (lux  equation  is  considered  where  the  waves  are  assumed  to  be  steady  state 
and  propagating  normal  to  the  shoreline  over  straight  and  parallel  depth  contours, 

3Sxx  ..    .-.    3n 


3x 


9(h  +n)    8^  (5.1) 


Equation  5.1  indicates  that  the  change  of  the  excess  momentum  flux  due  to  wave 
motion  in  the  cross-shore  direction  is  mainly  balanced  by  the  mean  water  slope  dfj/dx, 
assuming  no  gradients  in  the  Y-dircction  and  neglecting  the  frictional  forces. 

It  is  common  to  separate  the  nearshore  zone  into  two  dynamical  regions  of 
"outside"  and  "inside"  the  surfzone. 

(1)  Outside  the  Surf  Zone  : 

If  waves  are  assumed  to  propagate  onshore  without  any  energy  loss,  then  the 
onshore  component  of  the  radiation  stress  is  expected  to  increase  steadily  as  the  depth 
decreases.  Therefore,  the  mean  water  level  is  lowered  (set  down)  by  the  presence  of  the 
unbroken  waves.  Longuet-IIiggins  and  Stewart  (1962)  integrated  equation  5.1  by 
assuming  that  rf  <  <  h  to  obtain 


—  1  a2k 

n      -      "    2    sinh    2kh 


(5.2) 
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This  equation  is  second  order  in  local  wave  amplitude.  The  setdown  increases  steadily 
until  the  breaking  zone.  Without  referring  to  the  radiation  stress,  Longuct-lliggins 
(1967)  also  derived  the  above  solution  using  Bernoulli's  integral  equation.  This  simple 
technique  will  be  illustrated  and  applied  below. 

(2)  Inside  the  Surf  Zone 

As  waves  break,  the  wave  amplitude  inside  the  surf  zone  is  assumed  to  be 
controlled  by  the  total  local  water  depth  including  the  wave  setup  rj  ,  simulating  a 
spilling  breaker  type,  i.e. 


=      y  (h  +n)  (5.3) 


The  cross-shore  component  of  the  radiation  stress,  Sxx  ,  is  a  function  of  the  energy 
density  which  decreases  inside  the  surf  zone  and  reduces  to  3/2  li  in  shallow  waters 
such  that 


3         2  -x2 

S  =      TP9Y     <h  +n)  (5.4) 

xx  4 


Both  the  gradients  of  SYY  and  the  forcing  in  equation  5.1   have  a  negative  slope. 

\  A 

Therefore,  the  change  in  the  MWL  inside  the  surf  zone  is  expected  to  be  positive, 
representing  a  wave  setup.    Integrating  equation  5.1  gives 


=       N(hb   -h)     +    r,b  (5.5) 


7        1 

where  N  =  (1+   2/(3  y    ))'    .    The  subscript  b  indicates  the  value  at  the  breakerhne 

where  ^  =•  y  JL  /8  is  determined  from  equation  5.2  .  Equation  5.5  shows  an 
increase  in  the  mean  water  level  with  decreasing  depth,  balancing  the  decreasing 
radiation  stress  due  to  wave  dissipation  inside  the  surf  zone.    It  should  be  noted  that 


Bernoulli's  integral  outside  the  surf  zone  is  not  applicable  inside  the  surf  zone  since  the 
motion  is  no  longer  irrotational. 

Bovven  et  al.  (1968)  demonstrated  the  above  theoretical  results  by  careful 
measurements  in  the  laboratory.  Experimental  results  agreed  well  with  the  theory. 
However,  just  outside  the  breaking  point,  their  experimental  values  were  less  than  the 
theoretical  predicted  since  the  theory  does  not  consider  energy  dissipation.  Also,  very 
close  to  the  shoreline,  Bowcn  et  al.  (1968)  found  a  residual  wave  height  (possibly  due 
to  a  standing  wave)  to  persist  and  corrected  the  above  equation  by  letting  the  setup 
slope  be  tangential  to  the  beach  slope,  i.e.  dr|/dx  ->  -  dh/dx  as  (h  +  r\)-^  0  .  Van 
Dorn  (1976)  found  that  the  setup  gradient  at  the  shoreline  is  actually  steeper  than 
predicted  by  equation  5.5,  which  agrees  with  Bowen  et  al.  (196S).  Battjcs  (1974)  and 
Battjcs  and  Jansscn  (1978)  used  a  random  wave  model,  which  predicts  smaller  setup 
values. 

In  addition  to  the  wave  setup  laboratory  investigations,  a  few  field 
measurements  have  been  reported.  Dorrestcin  (1961)  measured  the  change  in  mean 
water  level  across  the  surf  zone  during  times  when  the  significant  offshore  wave  height 
ranged  from  0.8  to  1.6  m.  lie  measured  a  maximum  setup  of  0.15  m.  Guza  and 
Thornton  (1981)  measured  the  setup  at  the  shoreline  when  the  significant  offshore 
wave  height  varied  from  0.6-1.6  m.  The  maximum  setup  at  the  shoreline  was  found  to 
be  about  0.17  1  Iq  ,  where  11^  is  the  deep  water  significant  wave  height.  Ilolman  and 
Sallcnger  (1985)  measured  the  setup  under  conditions  of  incident  significant  wave 
height  varying  from  0.4-  4.0  m.  They  found  the  setup  to  vary  linearly  with  the  surf 
similarity  parameter  ^()  =  P  (Hq  /Lq  )    '     . 

2.  Changes  in  Setup  and  Set  down  Due  to  the  Presence  of  Long  Waves 

Changes  in  the  wave  setup  and  setdown  due  to  the  presence  of  long  waves  can 
be  demonstrated  by  considering  the  same  wave  field  used  before,  i.e.  oblique  short 
waves  riding  on  a  normally  incident  long  period  standing  wave.  The  wave  setdown  is 
evaluated  applying  Bernoulli's  integral  equation  (Longuet-IIiggins,  1967).  first,  the 
vertical  momentum  equation  is  integrated  over  depth  and  averaged  over  time. 


(p  +pw2)  |       -    pgn      =      0  (5.6) 
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where  p  is  the  pressure  and  the  ovcrbar  denotes  mean  quantities.  Next,  the  time 
averaged  Bernoulli's  integral  equation  is  evaluated  at  z  =  0,  assuming  the  flow  is 
irrotational, 


1     rTT2    4.    w2    +   w2)  I  +    Constant    =    0 


P|  +   ±p(Ir    +   v      +  w")  |  +   Constant   =    u  (5-7) 

z  =  0 


2  z  =  0 


The  pressure  term  is  eliminated  by  combining  equations  5.6  and  5.7  to  obtain 


-     =      -    -L(u2    +   v2    -   w2)  |  +   Constant  (5  g) 


The  difference  in  mean  sea  level  5  \\  is  obtained  by  applying  equation  5.8  at  two 
different  locations  (xj  ,yj  ,0)  and  (X2  ,}'2  .0).  which  avoids  evaluating  the  constant. 


1 

-     =      _      l(u2    +   v2    -   w2)|  I  (5.9) 

on  2g  z=0      2 


For  the  combined  short  and  long  standing  waves,  equation  5.9  is  used  to  evaluate  the 
sctdown  as  follows,  liach  velocity  component  is  expressed  as  the  sum  of  the  short  and 
long  standing  wave  velocities  as  given  before  in  chapter  2.  Let  location  2  be  in  deep 
water  (h  -^  °°  )  where  the  sctdown  is  zero.  Substituting  into  equation  5.9  and 
averaging  over  the  long  wave  period  gives 


a2    k  —        — 

-     _  s    K  1,2  2.  (5.10) 

n      "      "    2    sinh    2kh   "    2?(u£    "   V 


where  subscripts  s  and  1  refer  to  short  and  long  waves.  The  first  term  in  equation  5.10 
is  identified  as  the  the  sctdown  solution  (equation  5.2)  when  considering  the  short 
wave  only.  At  the'  brcakerlinc,  the  depth  controlled  breaking  model  is  utilized  and 
sctdown  reaches  a  maximum  to  give 
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2  2 

Y    (h.    +  n0,  ) 


b        'g,b'  1,22.  ,-  f.„ 

nb nr 2^(u^  -  V  (5J1) 

b 


where  the  subscript  b  refers  to  conditions  at  the  breakerline. 

Wave  setup  can  be  evaluated  by  spatially  integrating  the  X-momentum  flux 
equation,  where  the  radiation  stress  component  Sxx  for  the  combined  wave  fields  is 
given  by  equation  3.74  .  The  constant  of  integration  is  evaluated  by  applying  the 
obtained  setup  equation  at  the  breakerline  where  the  change  of  the  mean  sea  level  is 
known  from  equation  5.11  .  The  resulting  equation  is  quadratic  in  i|  and  has  two 
roots.  One  root  is  disregarded  since  it  unreasonably  predicts  setdown  at  the  shoreline. 
The  second  root  gives 


=      -(h   + 


A2   -Q    )  (5-12) 


where 


1,,    2/u2    _    u2 

b    "'    "  wb "b 


-In*    +    2hh-h    *i{*     (hb-h    ) 


+   a2(V  "    1)  (Jn    (X)    '   J0(X)) 


(5.13) 


i 


b 


+    2a2(jl    (X)     -    jJ(X))}] 


Using  Longuct-lliggins'  setup  model  for  short  waves  as  a  reference,  the  setup  is 
relatively  increased  closer  to  the  shoreline  due  to  the  presence  of  the  long  standing 
wave  (fig  5.1  ).  However,  the  setdown  is  also  increased  at  the  breakerline  compared 
with  the  reference  model.  The  relative  increase  in  setup  appears  to  be  over-predicted 
and  suggests  the  importance  of  including  damping  terms,  such  as  frictional  dissipation 
and/or  percolation,  in  the  swash  zone. 
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Figure  5.1     Wave  setup  due  to  both  incident  and  infragravity  waves 
compared  with  Longuet-Higgins  and  Stewart  model  (1962). 
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B.       LONGSHORE  CURRENT 
1.  Background 

When  waves  obliquely  approach  the  surf  zone,  a  mean  current  parallel  to  the 

shoreline  is  generated,  which  is  important  in  transporting  sediment  alongshore.    In  the 

simplest  formulations,  the  longshore  current  is  assumed  steady  and  two  dimensional, 

being  independent  of  the  longshore  topographic  changes.   The  linear  Y-momentum  flux 

equation  is  simplified  such  that  the  gradient  of  the  radiation  stress  component  d  SVY  /dx 

yx 

is  balanced  by  the  bottom  shear  stress, 


Sjx     =      _  (5.14) 


ox  y 

Outside  the  surf  zone,  the  wave  energy  is  assumed  non-dissipative  and  the  excess 
momentum  flux  in  the  Y-dircction  is  a  conservative  quantity.  Therefore,  no  driving 
force  is  expected  outside  the  surf  zone  for  generating  a  longshore  current  (Bowen, 
1909).  Inside  the  surf  zone,  energy  is  dissipated  due  to  wave  breaking,  causing  a 
change  in  the  momentum  (lux.  The  onshore  changes  in  momentum  llux  drive  a 
longshore  current  primarily  confined  to  the  surf  zone. 

Numerous  models  of  longshore  current  on  a  long  straight  beach  have  been 
done  since  the  introduction  of  the  radiation  stress  concept.  For  obliquely  incident 
waves  on  plane  beaches,  analytical  solutions  were  derived  for  monochromatic  waves  by 
Bowen  (1969)  and  Longuet-I  Iiggins  (1970  a,b)  and  for  random  waves  by  Thornton  and 
Guza  (1986).  Numerical  treatments  (e.g.  Thornton,  1970;  Jonsson  et  al.,  1974; 
Madsen  et  al.,  1978);  Wu  and  Liu,  1984,  and  others)  have  employed  arbitrary  beach 
profile  and/or  including  setup  and/or  non-linear  inertial  terms.  In  monochromatic 
wave  models,  lateral  mixing  is  introduced  to  smooth  out  the  discontinuity  in  the 
longshore  distribution  due  to  the  intensive  current  shear  at  the  breakcrline  (Bowen, 
1969;  Thornton,  1970;  Longuet-I  Iiggins,  1970a).  For  random  waves  models,  there  is 
no  need  to  include  the  lateral  mixing  (Collins,  1970;  Battjes,  1972)  since  the  waves 
transform  and  break  over  a  considerable  horizontal  distance  resulting  in  a  smooth 
change  in  the  rate  of  energy  dissipation  and  S  x  . 

Periodic  variation  and  unsteadiness  in  the  longshore  current  measurements 
have  been  reported 'by  Wood  and  Meadows  (1975)  and  Meadows  (1976),  but  neither 
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analytical  nor  numerical  solutions  have  been  attempted.  Meadows  (1976)  suggested 
three  distinct  longshore  velocity  components,  a  steady  component  and  two  fluctuating 
components,  with  one  at  a  short  wave  frequency  and  the  other  at  low  frequency. 
Guza  and  Thornton  (1978)  and  Holman  and  Bowen  (1984)  observed  low  frequency 
components  in  the  longshore  velocity  spectra.  Since  the  longshore  current  is  a 
timc-avcraged  velocity,  Guza  and  Thornton  (1978)  stated  that  an  appropriate  temporal 
averaging  time  for  mean  longshore  current  is  unknown.  The  above  studies  suggest  that 
the  low  frequency  wave  motion  can  be  important  inside  the  surf  zone. 

Infragravity  waves  are  included  in  the  following  wave  description  inside  the 
surf  zone,  allowing  a  more  complete  description  of  the  wave  field.  The  changes  in  S 
due  to  obliqueness  of  the  short  wave  and  the  total  bed  shear  stress  are  derived.    Then, 
the   steady   and   unsteady   longshore  current   formulations  are   considered   on   plane 
sloping  beaches. 

2.  Wave  Refraction  and  Radiation  Stress 

The  depth  contours  are  assumed  to  be  straight  and  parallel,  and  obliquely 
incident  waves  are  therefore  refracted  according  to  Sncll's  law, 


sin    ct 

sin    a      =      £•     =      Constant  (5.15) 

c  c 


where  c  is  the  phase  speed.    The  S..Y  component  for  the  combined  short  and  long 

yx 

waves  is  given  by  equation  3.77  .  Assuming  a  small  angle  of  wave  incidence  (cos  a  ~ 
1),  Sncll's  constant  (sin  a  /c),  which  contains  all  the  wave  angle  information,  can  be 
introduced  into  equation  3.77  to  yield 


1       , s i n    a,      /— —      2  ,  u    ,     „    x  2 
yx 


S„„      =      2pg(_ a- )     ^  Y     (h   +    V  (5.16) 


Another  approach  for  deriving  equation  5.16  can  be  shown  by  using  equation  3.67, 
where 
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EC 
c 

yx  C  (5.17) 


S  =      — 2-   sin    a   cos    a 


By  assuming  a  small  angle  of  wave  incidence,  5.17  is  reduced  to  E  c  (sin  a  /c)  where  E 
is  given  by  1/2  p  g  y2  (h+  i^j  )~  .  This  will  give  the  same  solution  as  equation  5.16  . 
This  result  can  be  attributed  to  the  fact  that  standing  waves  by  themselves  do  not 
contribute  directly  to  the  longshore  momentum  flux  since  they  are  oriented  in  the 
cross-shore  direction,  and  no  S  forcing  is  expected  for  a  normal  standing  wave.  The 
contribution  of  the  infragravity  waves  manifests  itself  by  the  generation  of  momentum 
at  the  side  bands.  Averaging  equation  5.16  over  the  long  wave  period  gives 


S    ..      =      t^9 


2 ,L2         1      2 T2 


yx 


2 


5-iH^)     /gh   Y    (h      +   |a£J0(X))  (5-18) 


That  is  to  say,  the  S  forcing  is  increased  by  the  formation  of  the  side  bands  (second 
term  in  bracket).  The  longshore  current  forcing  term  d  SVY  /  d  x  for  both  steady  and 
unsteady  cases  is  compared  in  Fig.  5.2  with  the  Longuct-I  liggins  (1970a)  model  where 
only  short  waves  are  considered.  In  general,  the  combined  forcing  by  short  and  long 
waves  increases  in  the  vicinity  of  the  shoreline  due  to  the  relative  growth  of  the  side 
bands. 

3.  Bottom  Shear  Stress 

Inside  the  surf  zone  the  bottom  stress  can  no  longer  be  neglected,  especially 
when  considering  the  change  of  momentum  flux  in  the  Y-direction.  Assuming  the 
quadratic  shear  stress  formula,  the  averaged  alongshore  bottom  shear  stress,  exerted  by 
the  oscillating  wave  motion  is  given  by 


t         =      p  cf  |uT|uT    3  (5.19) 


where  c^  is  the  friction  coefficient  and  u-p  is  the  vector  sum  of  the  current  and  the 
inviscid   orbital   velocities  for  both  short  and  long  waves  measured  just  above  the 
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Figure  5.2     Gradient  of  the  longshore  current  forcing  across  the 

surf  zone. 
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bottom  boundary  layer.  The  absolute  value  is  necessary  to  insure  that  the  stress 
reverses  with  velocity  reversals.  Longuct-Higgins  (1970a)  assumed  a  small  angle  of 
wave  incidence  and  a  weak  current  compared  to  the  water  particle  motion  of  the 
waves,  and  he  was  able  to  simplify  5.19  to 


t^       =      p  cf  |uT|    V  (5.20) 


where  V  is  the  mean  motion  in  the  longshore  direction  which  is  assumed  uniform  over 
depth.  The  weak,  current  assumption  has  been  shown  to  be  a  reasonable 
approximation  for  many  field  applications  (Liu  and  Dalrymple,  1978).  The  total 
velocity  vector  u-r  can  be  resolved  into  its  components  in  the  cross-shore  and 
alongshore  directions: 


uT      =       {(us  +uj     +    U)I    +    Uvs  +vj     +   V}J  (5.21) 


where  u  and  Ui  are  the  X-componcnt  of  the  orbital  velocities  and  i,j  arc  the  unit 
vectors  in  the  x  and  y  directions  respectively.  U  and  V  are  mean  velocities.  Assuming 
uniform  currents  alongshore,  the  mean  cross-shore  velocity  U  is  zero  according  to  the 
conservation  of  mass,  since  the  beach  forms  a  boundary  in  the  X-direction. 
Considering  an  oblique  short  wave  riding  on  a  cross-shore  long  standing  wave,  i.e.  Vj 
=  0,  equation  5.21  can  be  rewritten 


u        =      {(u     cosa   +   u£)}i   +    {(us  sina)    +   vlj  (5-22) 


The  absolute  value  of  the  total  velocity  vector  is 


u    |       =       \/(u     cosa   +   u„)       +    (u     sina    +   V)  (5.23) 
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Assuming  a  small  angle  of  wave  incidence  (a  <  10°  ),  cos  a  can  be  approximated  to  be 
unity,  which  simplifies  5.23  to 


uj  \/(u      +    u    )2    +    V2    +    2u   V  sin    a  (5.24) 


It  follows  that  u    sin  a  is  small  and  can  be  neglected,  then 


/  (u      +    un)2    +   V2  (5-25) 


T1  V       s         «  2- 


which  can  be  rearranged  into 


u 


T 


~^-)2 


s    -    V      V1    "     'VV  (5.26) 


Note  that  us  +  Uj  =  Uj  i  according  to  equation  3.56  .   Then,  using  the  weak  current 
assumption,  V/u-r  <  <  1,  to  simplify  the  analysis  gives 

|uT|       =       |us    +    uj  (5.27) 


Since  the  two  velocities  us  and  Uj  are  almost  collincar  (both  acting  almost  in  the 
X-direction),  equation  5.27  gives  the  basic  relationship 


uTl      -      |usl    +    |ua 


(5.28) 


Equation  5.28  can  be  evaluated  using  3.57  and  2.17  and  averaged  over  the  short  wave 
period  to  give 


uTl     =    /57Ht%M  +  v0<x>l  +»1U.i«x)  «m  »»tj)         (sm) 
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The  alongshore  bottom  shear  stress  averaged  over  the  short  wave  period  is  then  given 
by  the  use  of  5.20  , 


T 


=      pcf    /g/hi^(|h+    a,J0(x)  |+a£  |J1(x)     sinwft|}V  (5.30) 


Since  5.29  is  still  time  dependent,  the  expression  is  averaged  over  the  long  wave  period 
to  fit  the  steady  formulation 


|uT,    =  2/71^/1^  {lYh  +  a1(yJ0(x)  +  J1(x)l}.  (5.31) 


The  corresponding  longshore  bottom  stress  is  given  as 


iy    =    |pcf  •57E{|7h+  a^trVxH+IVx)  D>v  (532) 


4.  Longshore  Current  Models 

Obliquely  incident  short  waves  riding  on  a  normally  incident  long  standing 
wave  arc  considered  to  develop  steady  and  unsteady  longshore  models. 

a.    Unsteady  Model 

Neglecting  the  nonlinear  mean  advective  terms,  the  linearized  Y-momcntum 
flux  equation  has  the  form 


)M  ciS, 


driv         "~YX  -    i  (5.33) 


where  M     =  V  p  h  is  the  mean  longshore  mass  transport  averaged  over  the  short  wave 
period.    No  variation  in  the  longshore  direction  is  allowed  in  equation  5.33,  and  the 
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lateral  mixing  is  ignored  for  simplicity.  Assuming  an  unsteady  longshore  current 
appears  to  be  more  general,  especially  with  the  presence  of  irregular  oscillatory  waves 
(Wood  and  Meadows,  1975).  Thus,  equation  5.33  can  be  written  in  terms  of  the 
longshore  current  V: 


3S 


p  h|v  +  p  Cf  |uT!v 


ax 


(5.34) 


The  onshore  divergence  of  the  radiation  stress  S      is  considered  the  significant  forcing 
for  generating  longshore  currents  and  is  given  by  differentiating  3.77  , 


t) 


jX 


(x,t) 


=      N*{B*(x)    +    CMxlcos    u)£t   +    DMxlcos    u^tJ     (5  35) 


where  N    is  a  constant  and  B  '  C    and  D    arc  functions  of  space  only.     N    is  given  as 


(5.36) 


•a     .sin    a.     , 
N*       =       p    -re- ( — )y 


j  tan  3    h 


5/2 


and 


=      a 


3/2. 


3h/gJQ(x)    -    2h"'  *"ojn  J,  (X)  - 


ri 


'gK 


D: 


r/gtanBJQ(X)    -    4  u^/hJ^X)" 


/gh 


The  longshore  current  forcing,  5  S      /<3  x  decreases  onshore,  then  increases  closer  to 

yx. 

the  shoreline,  having  a  singular  point  at  x  =   0  (see  Fig  5.2     );  the  computed  values 
stop  at  x  =   2  cm.    The  forcing  term  is  shown  to  oscillate  at  the  long  wave  frequency 
(Fig  5.2     ).    Dividing  equation  5.34  by  p  h  and  calling  the  resulting  forcing  term  q 
(x,t),  equation  5.34  can  be  simplified  to 


U      +    P*(x,t)V 

o  u 


=      q*(x,t) 


(5.37) 
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where  p    is  given  by 


p*(x,t)       =      c      vW^lh  +a£JQ(X)  I    +aJj1(X)    sin'u>£t|}     (5.38) 


It  can  be  shown  that  the  time  dependent  term  (second  term  between  brackets)  is  small 
compared  with  the  short  wave  contribution  in  the  swash  zone,  since  Uj  /  |ug  |  ->-  0  at 
the  shoreline.    liquation  5.38  then  reduces  to 


P.(x,tl       =      cf   v/g/h3{2Xih   +   ajlJ0(X)|]  (5-39) 

Equation  5.37  is  completed  by  specifying  the  initial  condition  V(x,0)  =  0.  An 
anal>tical  solution  to  the  unsteady  longshore  current  is  obtained  using  an  integrating 
factor.  Despite  the  fact  that  both  the  forcing  term  and  the  bed  shear  stress  have  a 
singular  point  at  the  shoreline,  the  resulting  longshore  current  is  well  behaved  and  has 
a  finite  value  at  the  shoreline.   The  general  solution  can  be  expressed  in  four  terms  as 


v(x'fc)       =      Vsteady    +   Vunsteady    1    +   Unsteady    2  (5.40) 

+    V 

transient 


The  first  term  is  a  steady  term  and  is  given  by 


2 
v  =       5tt  /Sin    n,     y      tan8      .2 

steady  8  (      c       '       c      A*      gh 

„2    2 

tt  .  s  i  n    a. .     '    a  l  ,  2  , 


(5.41) 


-    4oj   J    (x)  J    (x)     s^h] 
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where  A  =  y  (  h+  aj  Jq  (/))  .  The  steady  term  is  similar  to  the 
steady  solution  given  by  Longuet-Higgins  (1970a)  differing  only  by  the  second  term  in 
the  A  expression.  The  second  term  in  5.41  is  an  additional  steady  component  as  a 
result  of  the  side  band  formation.  The  two  unsteady  terms  arc  functions  of  the 
frequency  of  the  long  wave  and  its  first  harmonic  and  are  given  by 

'     T      9a?. sin    a} 

Unsteady    1      "      2  [4   ^    A*2    +    ^2  jy}         c 


2    3 

x     [3  tan6JQ(x)h    /gh   -    2u^J1(:A)h    ] 


x     (P*   cosoj.t    +    oo„    sinu.t) 

A.  JC  *> 


(5.42) 


2    2    2 

tt    v    a  „  g 

v  =    ! i_ (sin  a) 

unsteady    2  ._ r    2       ,.2  2      2.  2,         c 

32  [c      gA*       +    tt      u;,  h    ] 


[tan  3Jq(x)    ^h3    -    4    u)£JQ  (  X)  J^  (X)  h2] 


(5.43) 


(P*  cos2w    t    +    2u)?sin    2w    t) 

The  fourth  term  in  equation  5.40  is  the  constant  of  integration 
and  is  evaluated  by  applying  the  initial  condition  at  time  zero.  This  term  is  identified 
as  a  transient  term  which  grows  very  rapidly  to  a  steady  state  within  a  few  long  wave 
periods. 

)  I  le~P*t 

transient      =      "  [  (Vsteady  +  Unsteady    1  +Vunsteady    2      t=Q 


In  the  vicinity  of  the  shoreline,  the  longshore  current  is  found  to  be  finite  and  steady 
and  is  given  by 

tt         Y  tanB        q(sin    a) 
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The  spatial  distribution  of  the  longshore  current  is  shown  in  Fig  5.3  .  The  dominance 
of  the  side  bands  and  the  oscillating  components  starts  at  the  most  shoreward  1/5  of 
the  surf  zone  width.  At  the  breakerlinc,  intensive  shear  occurs  as  expected  in 
monochromatic  models.  Lateral  mixing,  if  added,  will  diffuse  the  shear  and  smooth  the 
longshore  distribution  as  described  before.  Comparison  with  Longuet-Higgins'  (1970a) 
model  shows  an  onshore  increase  of  the  longshore  current  and  a  relative  decrease 
closer  to  the  brcakerline.  The  decrease  is  due  to  the  the  increase  in  the  bottom  shear 
stress. 

The  time  variability  of  the  longshore  current  is  evident  in  Fig  5.4  ,  where 
the  current  is  shown  to  oscillate  at  the  infragravity  wave  frequency.  In  the  vicinity  of 
the  shoreline,  the  first  harmonic  (Vunsteadv2  )  dominates  the  fundamental  frequency 
component.  The  transient  term  is  just  an  initial  response  to  the  wave  system  that 
becomes  negligible  within  a  few  infragravity  wave  periods. 

b.   Steady  Model 

For  steady  flow,  equation  5.14  is  used  in  which  the  cross-shore  gradient  of 
Syx  is  balanced  by  the  bottom  stress.  As  in  the  analysis  of  the  unsteady  model,  wave 
setup  is  not  included.  Steady  solutions  with  and  without  considering  lateral  mixing  are 
developed.  Averaging  over  long  wave  periods  is  required  to  eliminate  the  time 
dependency  in  the  equations  describing  the  unsteady  models.  The  bed  shear  stress  is 
given  by  equation  5.32  ,  and  the  transverse  radiation  stress  component  is  obtained  by 
averaging  equation  5.35  over  the  long  wave  period.  Substituting  in  equation  5.14  ,  the 
steady  longshore  current  is  given  by 

2 

5tt    sin    q  >     y      tan  B  ,2 

steady  8  l      c       '       cfA**        9 


2    2 


-  4w£jQ  (;,)  jl(x)    /gh] 


(5.46) 


** 


where  A       =  y  h  +  aj  (y  Jq  (x  )  +  Jj  (X  ))•    The  first  term  of  5.46  is  similar  to  the 

** 

steady  term  in  the  unsteady  case  except  for  A      .    At  the  shoreline,  the  growing  side 
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Figure  5.3     Unsteady  longshore  current  distribution, 
at  time  (a)  t  =  long  wave  period  (b)  t  =   1/2  long  wave  period. 
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bands  will  have  a  steady  and  finite  contribution  to  the  longshore  current  given  by 
equation  5.45,  i.e.  for  both  steady  and  unsteady  models  the  magnitude  of  the  current  is 
the  same  at  the  shoreline. 

Ignoring  the  exchange  of  momentum  due  to  horizontal  turbulent  eddy 
transport  increases  the  current  shear  at  the  brcakerline  and  results  in  a  sharp 
discontinuity  in  the  velocity  distribution.  The  presence  of  any  horizontal  mixing 
couples  the  adjacent  water  elements  together  and  diffuses  the  momentum  in  the 
cross-shore  direction.  Adopting  the  concept  of  eddy  viscosity,  Longuet-lliggins 
(1970b)  expressed  the  mixing  term  to  be  proportional  to  the  mean  current  V  such  that 


3     ,         3V.  (5-47) 

T  =       - — (v    D— ) 


VX  n  o-A 


where  v  =  Np  x  J  gD  is  the  horizontal  eddy  viscosity,  D  is  the  total  water  depth  and 
N  is  non-dimensional  coefficient,  (0  <  N  <  0.16).  Thus,  the  steady  Y-momentum  flux 
equation  becomes 


Y2S      =       -    r       +    —(v   Dii-) 


V>  (5.48) 


3x  Y         ox  dx 


Equation  5.48  can  be  rewritten  as  a  second  order  differential  equation  in  V  such  that 


d  S      /  d  x  =  mj  (x)  d2   V/  d  x2  +  m2  (x)  d  V/  d  x  +  m3  (x)  V  (5.49) 


where 


mj  =  N  x   h  V  gh 


m2  -  5/2  N  h  Vgh 
m3  =  -  p  cf  |  u^  | 
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Figure  5.4    Time  variability  of  the  longshore  current 
(a)x  =  2.0m(b)x  =   1.0  m. 
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and  |uj  |  is  defined  by  equation  5.31  .  Equation  5.49  is  solved  numerically  using  a 
finite  difference  scheme.  The  coefficients  m}  ,  m2  and  m3  are  evaluated  at  mid-points. 
A  tri-diagonal  matrix  is  generated  over  a  domain  equal  to  twice  the  surf  zone  width. 
Two  boundary  conditions  are  specified  at  the  shoreline  and  away  from  the  breakerline. 
At  the  shoreline,  the  eddy  viscosity  term  goes  to  zero,  i.e.  the  mixing  length  goes  to 
zero  at  a  boundary,  and  equation  5.48  reduces  to  5.14  .  Therefore,  the  steady  solution 
without  considering  eddy  viscosity,  equation  5.45,  is  used  to  specify  the  boundary 
condition  at  the  shoreline.  This  approach  is  necessary  because  equation  5.49  is 
singular  at  x  =  0,  and  care  must  be  exerted  during  the  numerical  formulation  to  avoid 
mathematical  inconsistency.  The  longshore  current  is  assumed  to  vanish  at  a  distance 
of  one  surf  zone  width  seaward  the  breakerline. 

Steady  longshore  current  models  induced  by  short  and  long  waves  are 
shown  in  Fig.  5.5  compared  with  Longuet-Higgins  (1970a)  model.  A  comparison  with 
Longuet-Higgins  (1970a)  model  shows  a  relative  decrease  in  the  current  seaward  the 
inner  surf  zone.  This  can  be  attributed  to  the  increase  in  the  total  horizontal  velocity 
associated  with  the  bottom  shear  stress.  It  can  be  seen  that  including  the  lateral 
mixing  smooths  the  longshore  current  profile  at  the  breakerline  and  shifts  the 
maximum  velocity  onshore.  The  longshore  current  distribution  is  similar  to  the  results 
obtained  from  the  EOF  analysis  by  Guza  et  al.  (1986)  showing  a  non- vanishing  current 
at  the  shoreline  (Fig.  1.4). 

C.       LONGSHORE  SEDIMENT  TRANSPORT 

Incident  waves  arriving  obliquely  to  the  breaking  zone,  releasing  their  energy  and 
momentum,  generate  a  longshore  current  which  inturn  drives  the  longshore  sediment 
transport  along  the  beach.  Knowledge  of  sediment  transport  rate  is  essential  to  the 
understanding  of  various  coastal  engineering  problems,  such  as  topographic  changes  in 
beach  profiles  and  the  associated  marine  protective  works.  For  estimating  the 
longshore  sediment  transport,  there  are  three  general  approaches:  the  wave  power 
model,  the  traction  approach  and  Bagnold's  energetics  model.  In  the  wave  power 
model,  the  rate  of  total  longshore  transport  Q^  is  simply  assumed  to  be  proportional  to 
the  longshore  component  of  wave  powrer  ?o  at  the  breakerline  to  some  power  n  (Watts, 
1953), 
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Figure  5.5     Steady  longshore  Current  models  ,  with  and  without  adding 
lateral  mixing,  compared  with  Longuet-Higgins  model  (1970a). 
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Q*      •      K   P"  (5.50) 


where  K  and  n  arc  empirical  coefficients.   P*  is  given  by 


p;.      =       (ECg}b  COSab    sin   ab  {5'5° 


which  can  be  expressed  in  terms  of  the  radiation  stress  (Longuct-Higgins,  1972) 


P,      -      Syx   C  (5.52) 


Inman  and  Bagnold  (1963)  found  n  in  equation  5.50  to  be  unity  and  suggested  the  use 
of  the  immersed  weight  sand  transport  rate 


Iz      =       (Ps    "    Pw)9  a'Q£  (5.53) 


where  ps  ,  pw  are  the  density  of  sand  and  water  and  a'  is  the  sand  porosity  (roughly 
0.60).  Komar  and  lnman  (1970)  measured  longshore  sand  transport  rate,  combined 
their  data  with  other  data  to  seek  a  general  formula.  Most  of  the  data  falls  onto  a 
single  trend  satisfying  the  simple  relation 


I         =       0.77    P£  (5-54) 


-l 
where  P^  is  in  watts/m  when  specifying  wave  height  in  terms  of  Hrms  and  Qgis  in  m 

/day.    However,  liquation  5.54  is  empirical  with  no  real  consideration  to  sand  transport 

mechanisms.    In  addition,  the  coefficient  .77  has  a  high  degree  of  uncertainty. 


109 


The  traction  approach  includes  work  done  by  Bijker  (1971),  Swart  (1976), 
Madscn  and  Grant  (1976)  and  others.  For  instance,  Madsen  and  Grant  (1976)  adapted 
the  Hinstein-Brown  sediment  transport  equations  to  time-varying  sand  movements 
under  combined  waves  and  unidirectional  currents.  The  bottom  stress  and  resulting 
sand  transport  vary  with  time.  Their  model  is  relativly  complicated  and  requires 
parameters  for  which  no  data  are  available. 

Bagnold's  energetic  model  assumes  that  the  total  sediment  transport  is  the  sum  of 
the  bed  load  and  the  suspended  load  transport.  Bed  load  is  maintained  by 
grain-to-grain  contact  within  10  cm  of  the  bottom  while  the  suspended  mode, 
containing  generally  finer  grains,  is  lifted  from  the  bed  by  the  turbulent  action  of  the 
wave-induced  vertical  motion.  The  suspended  sediment  grains  are  supported  via 
turbulent  diffusion.  They  move  with  nearly  the  local  water  velocity.  However,  they 
are  also  falling  vertically  relative  to  the  local  supporting  fluid  with  a  fall  velocity,  W. 
Thus,  part  of  the  available  energy  is  extracted  to  maintain  the  particles  in  the  water 
column.  Bagnold  (1963)  developed  the  total  load  sediment  equation  for  steady, 
two-dimensional  stream  How, 


l      = 


lB    +    lS 

J  tan    <p    -    tan    3         w 

u 

-    tan    3  1 

(5.55) 


where 

i  =  total  immersed  weight  sediment  transport  (bed  and  suspended) 

Cjj  =  bed  load  efficiency 

Co  =  suspended  load  efficiency 

(j)  =  internal  angle  of  repose  of  the  sediment 

tan  (5  =  slope  of  the  stream  bed 

u  =  mean  velocity  of  the  stream 

W  =  fall  velocity  of  the  sediment 

il  =  rate  of  energy  production  of  the  stream 
Inside  the  surf  zone  (oscillatory  flow),  Bagnold  (1963)  relates  the  rate  of  immersed 
weight  transport  to  the  work  done  by  waves  and  current  on  the  bed  surface.    The 
mechanics  of  sand  particle  movement  are  described  as  a  back-and-forth  motion  by  the 
oscillatory  wave  motion  close  to  the  bed  with  essentially  no  net  transport.    Wave 
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energy  is  expended  in  supporting  and  suspending  the  moving  sand  above  the  bed. 
Once  the  sediments  are  suspended,  the  presence  of  any  unidirectional  current 
superimposed  on  the  to-and-fro  motion  can  then  produce  a  net  drift  of  sediment. 
Bagnold  (1963)  derived  the  following  relationship 


^6  (5.56) 


m 


where  \q  is  the  time  averaged  immersed  weight  transport  rate  per  unit  width  in  the 
direction  0  determined  by  the  unidirectional  current  uq  ,  H  here  is  the  available  wave 
power  due  to  energy  losses  by  the  fluid  drag  at  the  bed  surface,  um  is  the  magnitude  of 
the  orbital  velocity  of  wave  motion,  and  K'  is  a  dimensionless  coefficient.  The  success 
of  the  wave  power  and  Bagnold's  energetics  approaches  is  due  to  their  relative 
simplicity,  requiring  only  one  or  two  parameters  specified  from  the  data. 

Equation  5.56  has  been  used  as  a  basis  for  the  development  of  a  number  of 
longshore  transport  models.  Inman  and  Bagnold  (1963)  specified  the  various 
parameters  of  waves  at  the  brcakcrline  to  calculate  the  longshore  immersed  weight 
transport  rate, 


vb 

I        =      K'(EC    hcos    a      —  (5.57) 

*■  ^  m 


Equation  5.57  can  be  viewed  as  a  general  relationship  because  it  does  not  specify  the 
forcing  for  the  longshore  current.  It  could  be  due  to  a  tidal  current,  a  nearshore  cell 
circulation,  wind  or  to  oblique  wave  approach.  Komar  and  Inman  (1970)  required  K' 
to  be  equal  to  0.28  to  agree  with  their  data. 

Bowen  (1980)  applied  Bagnold's  model  to  the  problem  of  on-offshore  sediment 
transport  on  beaches  with  normally  incident  waves  and  onshore  steady  current  (no 
longshore  currents).  Bowen  used  Stokes  second  order  wave  theory  to  predict  the 
equilibrium  beach  profile  as  a  function  of  incident  wave  characteristics.  Using  a  similar 
approach,  Bailard  and  Inman  (1981)  derived  a  general  bed  load  equation  in  the 
longshore  and  cross-shore  directions  by  considering  time-varying  flow  over  a  sloping 
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bottom.  Their  solution  reduces  to  equation  5.56  for  steady,  two-dimensional  flow  and 
assuming  a  weak  current  and  a  small  angle  of  incidence.  Bailard  (1981)  extended  the 
work  of  Bailard  and  Inman  (1981)  to  include  the  suspended  sediment  load  to  give  a 
more  general  expression.  The  total  longshore  sediment  transport  of  Bailard  (1981) 
solution  is  given  by 

(5.58) 


+  pc,  -4  hi3  v 

where  the  first  term  on  the  R.H.S.  represents  the  bed  load  and  the  second  term  is  the 
suspended  load.  The  total  average  sediment  transport  can  be  obtained  by  integrating 
equation  5.58  across  the  surf  zone  from  zero  to  the  breakerline  (x^  ). 


Xb 
I.      -        /        i„    dx  (5-59) 

0 


Calibration  of  the  model  can  be  done  using  the  wave  power  model. 

There  is  considerable  uncertainty  concerning  the  relative  importance  of  bed  load 
versus  suspended  load  within  the  surf  zone.  According  to  field  measurements, 
Komar(1978)  suggested  that  the  suspended  sediment  transport  rate  accounts  for  less 
than  20%  of  the  total  transport.  However,  Downing  (1984)  estimated  the  suspended 
load  on  a  dissipative  beach  to  account  for  as  much  as  45%  of  the  littoral  drift. 
Sternberg  et  al.  (1984)  measured  the  suspended  sediment  rate  at  Leadbctter  Beach, 
Santa  Barbara  and  found  that  the  measured  suspended  sediment  longshore  transport 
rate  is  equal  to  the  total  longshore  transport  rate  as  predicted  by  equation  5.57  .  This 
result  agrees  with  the  data  acquired  by  Kana  and  Ward  (1981)  during  a  storm. 
Sternberg  ct  al.  (1984)  found  maximum  suspended  sediment  occurcd  at  the  midsurf 
zone  region,  and  pointed  out  that  the  suspended  sediment  appears  to  be  high  in  the 
swash-backwash  region. 

Kraus  ct  al.  (1981)  measured  the  longshore  sand  transport  distribution  using 
multicolor    fluorescent    sand    tracer;    a    distinct    bimodal    distribution    was    obtained. 


112 


Maxima  in  transport  occured  in  the  swash  and/or  within  the  breaker  zone,  where 
suspended  load  transport  is  expected  to  be  relatively  high  due  to  the  turbulence.  White 
and  Inman  (1986)  measured  the  longshore  sediment  rate  in  Santa  Barbara  and  also 
found  a  bimodal  distribution  across  the  surf  zone  with  maxima  near  the  break  point 
and  in  the  swash  zone  (Fig  5.6  ).  They  estimated  the  longshore  suspended  load  to  be 
roughly  10-30%  of  the  total  longshore  transport  rate. 

The  above  results  show  a  large  volume  of  sediment  in  the  swash  zone  which  is 
not  predicted  by  any  available  model.  Waddell  (1973)  found  a  low-amplitude 
migration  of  bed  forms  which  migrated  consistently  downslope  suggesting  that  uprush 
sediment  transport  was  primarily  through  suspension  sediment.  Bradshaw  (1980) 
pointed  out  that  the  swash-backwash  collision  process  (low  frequency  waves)  is 
significant  and  results  in  temporary  suspension  of  a  large  volume  of  sediment  in  the 
roll  wave.  Wright  ct  al.  (1982)  demonstrated  the  dominance  of  the  suspended  sediment 
in  the  inner  and  swash  zone  using  a  turbidimeter  situated  10  cm  above  the  bed.  The 
time  series  shows  pronounced  turbidity  maxima  at  infragravity  intervals.  In  general, 
the  predictions  of  littoral  drift  have  a  high  degree  of  uncertainty,  which  results  partly 
from  our  basic  inability  to  make  accurate  measurements  of  sand  transport  on  beaches. 

1.  Longshore  Sediment  Transport  due  to  short  and  long  waves 

The  incident  waves  arc  assumed  to  approach  obliquely  the  surf  zone  at  a  small 
angle  riding  on  a  cross-shore  long  standing  wave.  Following  Bailard  (1981),  equation 
5. 58  is  used  to  calculate  the  total  longshore  sediment  transport  distribution.  Based  on 
laboratory  and  field  data,  Bailard  (1981)  estimated  the  efficiency  coefficients  to  be  Cg 
=  0.13  and  e<j  =  0.025  and  W  =0.04  m/sec  .  The  total  horizontal  velocity  squared 
averaged  over  long  wave  period  is  used  in  evaluating  the  bed  load  term  in  equation 
5.58 

^      =      ^f2gh    +    a2(2){Y2j2(x)    +   ^^}  (5.60) 


Similarly,  the  averaged  cubic  total  horizontal  velocity  is  evaluated  to  calculate  the 
suspended  load  term, 
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uT|3       =      ^(273(gh)3/2    +    3Y2a£g/gh[2YJ0(X)    +J1<X)]  (5  61) 


2    2 


+    3v    a?    J0(x)g/g/h[J1(X)    +   2j    (  X)  ] 

+    a3(^)3/2[2f3J3(A)    +   J3(X)    +   ly2J2Q(X)J1(X)]} 

The  steady  longshore  current,  with  lateral  mixing  added,  is 
utilized  to  derive  the  sediment  alongshore.  The  longshore  sediment  transport 
distribution  due  to  the  combined  waves  is  shown  in  Fig  5.7  .  Note  that  including  the 
lateral  mixing  smooths  the  distributions  of  both  the  longshore  current  and  sediment 
transport  at  the  brcakerline.  The  resulting  distribution  qualitatively  simulates  the 
bimodal  distribution  given  by  White  and  Inman  (1986),  since  it  describes  the  maximum 
sediments  near  both  the  shoreline  and  the  brcakerline. 

Unlike  previous  models,  the  wave  energy  has  been  shown  to  exist  at  the  shore 
(in  the  swash  zone),  and  hence  observable  dynamical  processes  can  be  described 
analytically.  None  of  the  available  models  for  predicting  the  longshore  sediment 
transport  explain  the  increase  in  longshore  sediment  transport  in  the  swash  zone.  All 
the  previous  models  predict  no  transport  at  the  shoreline  which  has  been  shown  to  be 
unrealistic.  It  is  shown  that  the  longshore  sediment  transport  model  presented  here 
estimates  correctly  the  bimodal  distribution  and  shows  the  dominance  of  the  suspended 
load,  which  agrees  with  some  of  the  above  results.  The  main  improvements  of  the 
presented  model  are  the  inclusion  of  the  infragravity  waves  and,  therefore,  describing 
the  swash  dynamics  and  the  increase  of  sediment  transport  at  the  shoreline. 

D.       MODEL  VALIDATION 

The  wave-wave  interaction  approach  within  the  inner  surf  zone  suffers  from  our 
inability  to  model  the  breaking  process  accurately.  Nonetheless  it  can  lead  to  greater 
insight  into  the  inner  surf  zone  dynamics.  Several  assumptions  are  necessary  for 
simplifying  the  analysis  due  to  the  complexity  of  real  surf  zone  processes,  but  it  is 
believed  that  they  have  some  physical  basis.  It  is  assumed  that  the  linear  theory 
describes  the  waves  adequately  for  the  order  of  the  analysis,  and  inside  the  surf  zone 
spilling  breakers  are  assumed,  where  wave  height  is  strongly  controlled  by  the  local 
depth.  Non-breaking  long  period  waves  on  a  plane  sloping  bottom  are  assumed  to  be 
cither  leaky  mode  or  edge  waves  with  a  frequency  0)1  <  <  cos  .  Basic  assumptions  of 
straight  and  parallel  beach  contours  have  been  used  for  simplicity  where  no  longshore 
depth   variation   is  allowed.     Plane   sloping  bottom  is  assumed  for  simplicity.     The 
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Figure  5.6     Measured  longshore  sediment  transport  rate  i    distribution, 
Santa  Barbara  (alter  White  and  Inman,  1986/1 
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sensitivity  of  the  model  to  the  change  in  beach  slope  is  clear  since  infragravity  waves 
modify  the  local  water  depth  and,  therefor  modulate  the  total  energy.  Limitation  to 
narrow  banded  incident  and  infragravity  waves  is  required  for  the  model  input.  It  must 
be  noted  that  using  a  monochromatic  description  of  both  incident  and  infragravity 
waves  may  raise  questions  about  quantitative  predictions  since  the  processes  in  nature 
are  stochastic.  However,  the  developed  model  gives  a  number  of  qualitative  results 
which  may  be  useful  in  understanding  the  inner  surf  zone  dynamics.  Application  of  the 
model  to  different  nearshore  processes  leads  to  a  simple  physical  picture  of  the  nature 
of  the  wave-wave  interactions  in  the  inner  surf  zone.  The  validity  of  this  model  is 
confirmed  by  comparison  with  field  data,  where  the  resulting  distributions  at  least 
qualitatively  simulate  actual  distributions. 
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Figure  5.7     Longshore  sediment  transport  distribution  model  due  to  both 
short  and  long  waves  showing  a  bimodal  structure. 
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VI.  SUMMARY  AND  CONCLUSIONS 

The  modulation  of  incident  short  period  wind  waves  superposed  on  infragravity 
waves  are  investigated.  Modeling  the  short  waves  as  depth  limited  breakers,  the  short 
waves  are  shown  to  be  modulated  in  amplitude,  wavenumbcr  and  direction  due  to  the 
slowly  varying  depth  caused  by  the  long  standing  waves.  Furthermore,  modulation 
results  in  additional  energy  at  the  side  band  frequencies  near  the  incident  wave  peak. 
The  energy  calculations  show  a  gradual  growth  of  the  side  bands  onshore  with 
non-vanishing  energy  at  the  shoreline  (order  of  15%).  An  analytical  expression  is 
developed  to  describe  the  surface  elevation  of  the  modulated  short  waves  as  an  infinite 
sum  of  spectral  components  at  the  incident  and  side  bands  frequencies.  It  is  shown 
that  the  modulated  short  wave  transfers  energy  to  the  side  bands  as  the  depth 
decreases.  Therefore,  it  is  suggested  that  the  spectral  components  of  the  lower  side 
bands  may  strengthen  the  infragravity  wave  band  as  the  depth  decreases  resulting  in 
the  infragravity  waves  dominating  the  power  spectra. 

The  changes  in  the  wavenumbcr  and  direction  of  the  short  wave  due  to  the 
presence  of  a  long  wave  are  calculated  by  using  a  moving  frame  of  reference  that 
reduces  the  non-breaking  long  wave  to  a  slowly  varying  current.  On  a  plane  sloping 
beach,  the  wavenumbcr  modulation  is  shown  to  increase  with  increase  angle  of  wave 
approach.  The  excess  momentum  (lux  components  and  bed  shear  stress  due  to  the 
combined  waves  are  redcrived  analytically  and  show  periodicity  at  the  frequency  of  the 
long  waves. 

Spectral  analysis  shows  the  dominance  of  the  infragravity  waves  closer  to  the 
shoreline.  Using  narrow  band  data  with  high  resolution  (df  =  .004  Hz),  side  band 
peaks  can  be  identified  at  frequencies  (fs  ±  f«  ).  On  almost  all  days  and  for  all  the 
instruments,  the  upper  side  band  is  well  defined  in  the  power  spectra,  but  the  lower  side 
band  is  suppressed.  The  reason  for  the  lack  of  the  lower  side  bands  on  most  days  is 
not  known.  However,  breaking  of  the  low  frequency  waves  and/or  radiation  of  energy 
by  outgoing  component  of  the  standing  wave  may  cause  suppression  of  the  lower  side 
bands. 

Most  days  show  a  positive  correlation  between  the  incident  wave  envelope  and 
low  frequency  motion  inside  the  surf  zone,  while  a  negative  correlation  is  observed  in 


118 


deeper  water  in  agreement  with  Longuet-Higgms  and  Stewart's  (1962)  explanation  of 
forced  waves.  The  positive  correlation  is  explained  by  the  short  wave  modulation  due 
to  the  presence  of  long  waves.  However,  the  values  of  the  measured  cross  correlation 
coefficients  are  found  to  be  low  (although  statistically  significant),  partly  due  to  the 
reflected  wave -field,  that  decreases  the  coherence  (Elgar,  1985).  On  some  days,  the 
observed  zero  cross-correlation  was  found  at  approximately  the  breakerline  and  other 
days  in  the  inner  surf  zone.  This  may  suggest  a  nodal  point  for  a  standing  wave  and,  or 
a  transition  location  between  forced  and  free  long  waves. 

It  is  hypothesized  that  some  of  the  often  observed  discrepancies  with  present 
short  wave  dynamical  models  occurring  in  the  inner  surf  zone  are  due  to  the  interaction 
between  the  incident  waves  and  the  infragravity  waves.  The  derived  expressions  for 
energy  and  momentum  fluxes  are  utilized  to  study  wave  setup,  longshore  current  and 
sediment  transport  in  the  inner  surf  zone. 

Wave  setup  calculations  show  an  increase  in  the  mean  water  level  compared 
with  setup  due  to  the  short  wave  only.  This  is  because  the  breaking  short  waves  are 
depth  controlled  and  the  depth  is  slowly  varying  due  to  the  presence  o[  the  long  waves. 
Consequently,  energy  at  side  bands  is  generated  and  the  setup  is  increased. 

A  steady  longshore  current  driven  by  oblique  short  waves  riding  on  cross-shore 
long  period  waves,  is  found  to  be  finite  at  the  shoreline  due  to  the  presence  of  finite 
side  band  energy.  The  model,  in  general,  predicts  smaller  currents  than 
Longuet-Higgins'(  1970a)  model  because  of  the  relative  increase  in  bed  shear  stress  due 
to  the  presence  of  long  waves.  In  the  derived  analytical  expression,  contributions  by 
the  short  waves  and  side  bands  can  be  identified.  It  is  shown  that  the  cross-shore  long 
wave  does  not  contribute  to  the  forcing  term,  c  Svx    d  x  ,  since  it  is  averaged  out. 

It  is  more  general  to  consider  an  unsteady  longshore  current.  The  analytical 
solution  for  an  unsteady  longshore  current  shows  three  major  terms:  a  transient  term 
that  vanishes  after  a  few  long  wave  periods,  a  steady  term  similar  to  the  one  obtained 
in  the  steady  model,  and  an  unsteady  term  which  oscillates  at  the  long  wave  frequency, 
but  driven  by  the  momentum  changes  at  the  side  band  frequencies,  and  acts  closer  to 
the  shore  line.   This  unsteadiness  is  attributed  to  the  side  band  formation. 

An  improved  longshore  current  model  is  necessary  to  improve  the  longshore 
sediment  transport  formulation.  Based  on  the  sediment  transport  model  by  Bailard 
(1981)  and  using  the  developed  steady  longshore  current  model,  a  longshore  sediment 
transport  distribution  is  calculated.    The  model  predicts  sediment  transport  maxima  at 
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the  shoreline  and  the  breakerline.  The  results  Qualitatively,  the  results  compare 
favorably  with  various  field  data  by  Zenkovich  (1960);  Kraus  et  al.,  (1981),  and  White 
and  Inman(1986). 

It  is  recommended  to  include  the  wave  setup  in  the  longshore  current  models. 
Numerical  solution  for  the  unsteady  longshore  current  including  the  large  angle  of 
wave  incidence  and  assuming  a  strong  current  is  of  interest.  Numerical  solutions  can 
be  compared  with  the  obtained  analytical  results.  Actual  depth  profiles  are 
recommended  to  use  in  computing  the  longshore  sediment  transport  rate  to  compare 
quantitatively  with  the  data. 
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